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SUMMARY 


n. been developed which is capable of calculating the 

layer for the supersonic portion of a mixed- 
compression aircraft inlet operating at angle of attack. The calculatixin pro- 
cedure IS based on using a zonal solution algorithm. The supersonic core flow 
IS determined using the method of characteristics for steady three-dimensional 
ImI:. shock wave and the internal shock wave system are computed using 

procedure. The boundary layer flow adjacent to 
both the forebody/centerbody and the cowl is determined using an implicit finite 
difference algorithm. The flow in a shock wave-boundary le.er interaction re- 
gion IS, computed using an integral analysis. 

The computer program has the capability to compute the internal, flow field 

nntfnnM!’^h?“h Jt® 'discrete fitting of the internal shock wave system. The 
option in which the internal shock wave system is not fitted can be employed 
in situations in which the strength of the internal shock wave system is weak, 
continuities acceptable solution can be obtained by smearlBg the internal dis- 

<- 1 . ntol.ecular transport can be included in the computation of 

the f|bw about. the forebody, and in the computation of the internal 

flow in which the shock waves are not discretely fitted. This is accomplisJied 
J 4 .CX thermal diffusion terms in the governing partial 

differential equations as forcing functions, or correction terms, in the method 
of characteristics scheme. 

The boundary layer computation allows for the specification of the wall 
temperature or normal temperature gradient at the wall . --Moreover,, distributed 
wall b.leed effects can be accounted for. in the analysis. Laminar, transitional, 
or fully turbulent boundary layer flows Can be Calculated. 


The thermodynamic model, molecular transport properties, and turbulence 
model employed in the computer program are contained in a separate group. of 
subroutines. The assumed thermodynamic model is that-of. a thermally and calor- 
ically perfect gas. Dynamic viscosity is represented by Sutherland s law. 

Thermal conductivity is represented in terms of the dynamic viscosity and the 
laminar Prandtl number. Turbulent closure is achieved by use of a two layer 
eddy viscosity formulation based on mixing length and velocity defect concepts. ^ 
Alternct-ive formulations may be employed by modifying the existing subroutines 
or by replacing them. . 

The contours of. the centerbody and the cowl are represented by a separate 
subroutine. The existing subroutine has the capability to describe a variety 
of axi symmetric contours. Other • geometri es , such as those having elliptic or 
super-elliptic cross-sections, may be described by suitably_modifying the 
existing subroutine or by replacing it. _ 

A major assumption of the present analysis is that the cowl lip is contained 
in a given axial station. Moreov.er, it is assumed that both the centerbody 
contour and the cowl contour are smooth and have continuous first partial deriv- 
ati ves . 

The computer program cannot: 

1. compute a subsonic core flow, 

2. compute the external flow field about the forebody if the 

bow shock wave does not. exist entirely around the forebody, or 

3. compute the internal flow field if the bow shock wave is ingested 
into the annulus. , 

The computer program is written. in Fortran IV for. CDC 6000 Series, 7000 
Series, and CYBER computer systems. The program can be easily modified to be 
compatible with other computers. 
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SECTION I 
INTRODUCTION 


A computer program has been developed for calculating the flow field in- 
cluding boundary layer effects for the supersonic portion of a mixed-compression 
aircraft inlet operating at angle of attack. The general features of the. inlet 
geometry and the flow field are illustrated in Figure 1. The theoretical anal- 
ysis on which the computer program is based is presented by Vadyak and Hoffman 
(1). The present study represents an extension of an earlier investigation 
[(2), (3)] which was concerned solely with the computation of the supersonic 
core flow. 

This report presents a discussion of the computer program organization, 
descriptions of the subroutines, a discussion of the. input parameters, a brief 
interpretation of the output i.nformation, and sample cases to illustrate the 
application of the analysis. 


BOUNDARY LAYER 











SECTION II 

PROGRAM ORGANIZATION 


1 . INTRODUCTION 

In this section, the program overlay structure is presented, and some gene"- 
al comments concerning program input und parameter initialization are 
The subroutines are identified which contain the thermodynamic model ,^the t«^ahs- 
port property representations, the turbulence model, and the formulations used 
for representing the centerbody and the cowl contours. The flow field integra- 
tion-’options are briefly discussed, and spme. cpmme^^^^^^^^ program ouc- 

put are presented. 


2. INTEGRATION OPTIONS 

The computational flow regime is divided into two subregi 'es: the external 
flow field and the internal flow field (see Figure 1). If desired, only the 
external flow field may be computed.. Alternatively, if the solution is known 
at the axial station of the coXl lip. the internal flow field may be determined 
without employing the forebody flow integration option. 


The program has 
without the discrete 
core flow solution, 
might be employed if 
weak, and thereby an 
discontinuities. 


the capability to compute the internal flow field with or 
fitting of the internal shock wave system in the supersonic 
The option in which shock waves are not discretely fitted 
the strength of the internal shock wave system is relatively 
acceptable solution can be obtained by smearing all internal 


The program can be executed with or without invoking the three-dimensional 
boundary layer calculation. The boundary layer calculation allows for the 
specification of distributed wall bleed and prescribed wal^l temperature or 
norLl temperature gradient. Laminar, transitional, and fully turbient bound- 
ary layer flows may be computed. 

The analysis includes the influence of molecular transport in the compu- 
tation of the supersonic core flow by treating the viscous . 

Sion terms in the governing partial differential equations as ’ 

or correction termL in the method of characteristics scheme. At present^ the 
oroaram has the capability to include the influence of molecular transport in 
the^computation oAhe external flow field about the 

putation Of the internal flow field in which the shock waves are not discretely 
fitted. The program option In which the internal shock waves are discr ^ _ 
fitted does not have the capability to include the influence of "‘c^ccuUr trans- 
port in the computation, but rather assumes the flow to be inviscid.and adiabatic 
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3. PROGRAM jDVERLAY STRUCTURE 


The overall program consists of 121 routines (16 program routines and 105 
subroutines). The program Is too large to be stored continuously In computer 
memory; hence an overlay scheme Is employed. The overlay structure Is pre- 
sented In Figure 2. Three overlay levels are used: the resident overlay level, 
the primary overlay level, and the secondary overlay level. 

OVERLAY (0,0)* Is the resident overlay and controls the overall execution, 
of the program. OVERLAY (0,0) calls the primary level overlays^ (9VERLAY (1,0), 
OVERLAY (2,0), 0VERLAY (3,0), and (8VERLAY (4,0). 0VERLAY (1,0) is used for 
data input and parameter Initialization. jSVERLAY (1,0) calls the four second- 
ary level overlays, OVERLAY (1,1), 0VERLAY (1,2), 0VERLAY (1 ,3). and 0VERLAY . 
(1,4). 0VERLAY (1,1) Is used for the Internal generation of the suoersonlc 
core flow Initial data using an approximate Integration technique. '0VERLAY 

(1.2) Is employed for converting. the output from 0VERLAY (2,0) which generates 
the supersonic core flow Initial data using the Jones algorithm (4). 0VERLAY 
(^^3) is used to generate the three-dimensional boundary layer initial data for 
the forebody/centerbody geometry using the Adams finite difference algorithm (5). 
0VERLAY (1,4) is employed for -the Initialization of the boundary layer computa- 
tional parameters. 

The supersonic core flow solution Is computed using the primary, level over- 
lay, 0VERLAY (3,0), and the secondary level overlays, 0VERLAY (3,1 ), 0VERLAY 

(3.2) , and 0VERLAY (3,3). 0VERLAY (3,0) contains the preponderance of the unit 
processes and the Interpolation schemes. 0VERLAY (3,1) contains the control 
logic used In performing the Integration for the external flow field about the 
forebody. 0VERLAY (3,2) contains the control logic used In determining the 
Internal flow field in which shock waves are discretely fitted. 0VERLAY (3,3) 
contains the control logic used In determining the internal flow field In which 
shock waves are not discretely fitted. 

The three-dimensional boundary layer solution is computed using the primary 
level overlay, 0VERLAY (4,0), and the four secondary level overlays, 0VERLAY 
(4,1), 0VERLAY (4,2), 0VERLAY (4,3), and 0VERLAY (4,4). 0VERLAY (4,0) contains 
the overall control logic for the boundary layer computation. 0VERLAY (4,1) 
contains the routines used for the attachment line flow solution and the three- 
dimensional flow solution. 0VERLAY (4,2) Is used for the computation of the 
three-dimensional shock wave-boundary layer Interaction region. 0VERLAY (4,3) 

Is employed to account for boundary layer growth. 0VERLAY (4,4) initializes 
the boundary. layer solution on the cowl boundary using an approximate technique. 

The computer program has a restart capability In which an Internally gen- 
erated restart file Is retrieved from tape. When the restart option Is used, 
control Is. returned^ to the Integration control secondary level overlay which 
was in use at the time the Initial execution was terminated [e.g.. If the fore- 
body flow was being computed, control Is returned to 0VERLAY (3,1)]. 

★ 

Throughout the text, 0 denotes the letter 0. 
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4. DATA INPUT . 

The Ingut data are entered through both namelist input and fonnatted read 
statements. Eleven namelists must always be specified: namelists LISTl to 

LISTVl. Namelist LISTl contains the Input parameters which specify the_flow 
field Integration options, the number of planes of .symmetry, the output options, 
etc. Namelist LIST2. contains the parameters which .specify the.free-stream 
conditions, the .vehicle orientation, and parameters associated with the speci- 
fication of the Initial-value plane. Namelist LISTS specifies the number of 
circumferential and radial stations used In the computational mesh of-the super- 
sonic core flow solution. Namelist LIST4 '.specifies the thermodynamic model and 
the molecular transport properties. Namelist LISTS contains the parameters 
which specify the centerbody and. the cowl contours. Namelist LISTS contains 
the parameters which specify the boundary layer computational mesh. Namelist . 
LIST? contains the input parameters which specify the wall temperature or Its 
normal gradient for use In the boundary layer computation. Namelist LISTS 
contains the parameters used for the specification of distributed wall bleed. 
Namelist LIST9 Is used for specification of the turbulence model and transi- 
tion criteria. Namelist LISTIO Is used for specification of all convergence 
tolerances and Iteration limits. Namelist LISTll specifies the parameters 
used for debug output. 

In addition to the namelist. Input, the two-dimensional Initial-value plane 
for the supersonic core flow can be Input by formatted read statements. Alter- 
natively, If the forebody Is conical, ahead of the axial location. of the Initial- 
value plane, the Initial data may. be generated Internally In the computer program 
by one of two techniques. 

The Initial data for the boundary layer calculation can be Input by for- 
matted read statements. Alternatively, for Conical forebody flows, the boundary 
layer Initial data for the forebody/centerbody geometry can be generated Inter- 
nally in the computer program by a finite difference calculation. The cowl 
boundary layer Initial data can also be generated 1nterna.lly using an approxi- 
mate technlquBj^ 

All Input data are read In within routine LINKIO. After the data have 
been entered, selected Input parameters are checked for errors and consistency. 

If any errors are detected, appropriate messages are printed and the execution 
of the program -Is aborted. The input parameters are tested In subroutine INITIL. 

It should be noted that many Input parameters have default values and do 
not, have to be specified. The default values for most of the input parameters 
are spe_cif1ed in LINKOO. . 

After*all input parameters and data have been entered, parameter Ini tlaL- 
izat1on.(for most of the subroutines) Is performed In subroutine INITIL. 


5. INITIAL DATA SPECIFICATION 


The initial da.ta for the supersonic core flow must be specified on a space- 
like plane of constant x (see Figure V). The flow must be supersonic at every 
point specified. The x-coordinate axis constitutes the longitudinal axis of 
the centerbody and the. Cowl. Moreover! the mean. flow direction Is assumed to 
be In the x-coordlnate direction. 

If the forebody .flow field Is to be determined, the Initial -value plane 
must be specified at an axial (x) station that is upstream of the forebody compu- 
tational flow regime (see Figure.!). The solution is then found along the stream- 
lines that pass, through the data points specified on the Initial-value plane, 
although some streamline addition and deletion are performed on the ensuing so- 
lution planes. The axial location, of the last solution plane of the forebody 
flow field Is. adjusted to lie at the axial station of the cowl lip, and consti- 
tutes the Initial-value plane for the internal flow field computation. 

If only the Internal flow field Is to be determined, the Initial-value 
plane must be specified at the axial station which corresponds to the x-positlon 
of the cowl lip (see Figure 1). The cowl lip is assumed to be contained In a 
plane of constant x. For the Internal flow field Integration, a po.lnt redistri- 
bution is performed on the Initial-value plane. This point redistribution Is 
required to define streamlines which lie in the stream surface formed by the 
cowl boundary. The solution is then found along the streamlines that pass 
through the redistributed points on the plane at the cowl lip axial station. 

It should be noted that the internal flow field cannot be computed if the bow 
shock wave is Ingested into the annulus. . 

The supersonic core flow initial data may be specified externally by the 
user. The user supplied data are entered by a formatted read of the file I TAPI . 

Alternatively, if the forebody is conical ahead of the axial location of 
the Initial -val ue plane, the supersonic flow property field on the Initial-value 
plane can be generated internally in the computer program by one of two methods. 
One method of internally generating the initial-value plane is by using an ap- 
proximate technique which employs tne Taylor-Maccoll solution for the flow about 
a circular cone. at zero incidence. A superposition procedure. Is used to obtain . 
an approximation to. the flow. about a circular cone at nonzero angle of attack 
by neglecting the cross flow effects. This superposition procedure effectively 
amounts to computing the flow turning angle in the meridional plane of the. 
glven solution point, and then obtaining the flow properties at that point by 
applying the Taylor-Maccoll solution for a cone half-angle equal to the flow 
turning angle. It must be emphasized that this is only an approximate tech- 
nique, giving the well accepted Taylor-Maccoll solution at zero Incidence, but 
becoming Increasingly less accurate as the angle of attack Is Increased. Sub- 
routines LINKll, TMC0NE, INTGRT, L0GIC, C0NVRG, 0DESLV, and FUNCTM are used to 
genera.te Internally the supersonic core flow initialfvalue plane property field. 

The other method of Internally generating the supersonic flow Initial data 
Is based on the numerical Integration algorithm developed by Jones (4). The 
solution obtained by Jones for the flow about a Circular Cone at nonzero Inci- 
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dence has been well J°the^^n1tial-v 

conical up to the axial rtatlon where the initial h Jones algo 

Jones program Is the >^ecommended. source for (2.0). 

rlthm has been Incorporated Into the ^PMC, ZER0SN. FAl, FN, 

Routines LINK20-. DS0LV, and BODY comprise the Jones 

RKINIT. RKSTEP. owes, MAINM.MINIV. invert. ^ initial- 

algorithm and are f ®^.^; thf output^of the Jones algorithm compatible 

:?irtSi‘iicSSinu 

in subroutines LINK12* STENCL, REFLK, and J3 

If the boundary layer Jlow Is to J® fi®[^ is^to^be^ determined, 

be specified for the ®J"^ internal flow Is 

and for the forebody/centerbody flow and must be specified 

to be determined. Data ^^®,5®i4ii^^<n^aeneral not coincide with a plane 
on a body flr%Mc1fU of the torebody/center- 

of constant x. The axial ]°5t.i°?<iinc1dL with the axial location used for 

^St'sSrrsSSlo cfre Initial. Is%“lys l?“hrf1n? 

rpluSJJt? “staJL'^n’sirSS^in r l^Jce the cowl boundary layer thlcX- 
nesslls Identically zero at the cowl lip. 

The boundary layer /'»" '^‘ISLJtd’bra for^ntir«arrfTh"f1ll HM 
rinproft!;e=“fSfeifdyfcentl b 5y nltl^l d^^ and by a for«tted read of 
file ITAP3 for input of the .cowl. Initial data. 

Alternatively If the forebo^ Uc^^^^^^^^^^^ 

geneiated°1nterSally"l^ t^ ''(5f‘°uS”arSr'’tu?blllnt''tlOTs^^ 

sifsr^re 

;rTli:^S‘»Nrard%a ?»?r1se the Ad^ 

TiSm^U SUJ ^ 5^^1Ualu(Th^^p°pr^pr1ate p^ra.eters for the forebody/center 
body boundary layer computation. 

„3,„gTheTpir^XM5^^ 


6. THERMODYNAMIC MODEL 

The assumed thermodynamic iJP"*®;’ tncSrpoMted'lntrthrSmN^^ 

s:s s‘i.s f 

i>.*> 

modlflcatlohj^ 
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7. 


MOLECULAR TRANSPORT PROPERTIES 


funct4on^^f^h^temperat5?^i^^ assZed are assumed to be 

given by Sutherland's law 161 The lor-the dynamic 

the thermal conductivity is exDress»rf*^<n ii The assumed functional form for 
use Of the laminar Prandt nS^ viscosity thTouqh 

cosity and the thermal conductivity ffPcesentatlons for the dynamic vls- 
sultably modifying or reolaclna c^h. ;.'’® '"PPcporated -Into the analysis bv 
generation routines In 0VERLAy^l.3)\urre™?e ™3?fS«'on!''' 

8. TURBULENCE MODEL 

puter^Jogrlm^ o^dertfa^ f®®" Incorporated into the cora- 

on a mixing length formulation for the ^nJer^La^ J*’® "’°'‘®^ leased 

lation. for the outer reoion Th» mnwoi negpn and a velocity defect formu- 
for mass transfer at the wall Tho turbulence and account* 

pressed in terms of the eS5y viscositrLJ%"L%^'‘\\^®""®^ conductivitj is e?- 
ternate turbulence models may be incorooratAH^i^t^^!?!®"^ Prandtl number. Al- 
modifying or replacing subroitine TURb!^ ^ ^ analysis by suitably 


9. CENTERBODY AND COWL CONTOURS 


both the centerbody contour^and°the cowl^contou^"^ investigation assumes that 
purposes of geometry description thf a antisymmetric. For the 
of intervals, in any interval thA domain is divided into, a number 

lar input, or by supplying th^coeffi??entt^nlF^ ^P®t=ified by either tabu- 

function of X. It is assumed that hAth^+h^ a cubic polynomial written as a 
contour are smooth and have continumiQ ®®nterbody contour and the cowl 
etries, such as those havina ellifitiv derivatives. More arbitrary geom- 

into 

The initial-data aenlrat"^ replacing subroutines 
0VERLAY (1,3) and 0VERLAY ?2 ni ®“troutines in 0VERLAY 
tion for noncircular cross-sections *‘'''tP‘-AY <2,0) also require modifica- 


10. FLOW SYMMETRY 

gram ft?'’c«« ^HhlJ?; JhrbSmd^y lwSr’JSZ?!?'*'' P™- 

ntost general case is when no oianAc computation is not invoked. The 

compute the flew fl^d ?Sr ?ul f ?^3e.dlSvrL3?^“; ““d to 

second case is when one plane of SiJStrv «xi<?J ^ Incidence. The 

computing the flow field for axisyninetric inUt^At^f^^-^Pt^®" i"on 

case is when two planes of finw Inlets, at incidence. The third 

pute the -flow field for three-diSiSl^lJlh^ Iflh JP^^? “®®^ to corn- 

symmetry at zero- angle of attack The final ^wo planes of geometric 

synmetric. This option is uSd tn rALnl lu PJi°" i^ i^low is axi- 

Inlets at tero Incidmce’l 1" axlsynmetrlc 
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One flow symmetry option has been incorporated into the computer program 
for cases ia which the boundary layer is to be determined. This option is 
for the case of one plane of flow symmetry. 


11. OUTPUT 

Preliminary informc.tion is printed by LINKIO. The initial-value plane for 
the supersonic core flow and all subsequent solution planes are printed by sub- 
routine PRN0UT. The solution points on the space curves defined by the inter- 
section of the internal shock wave with the solid boundaries are output by 
subroutine SHKRFL. In addition to the position of^^and the dependent variables 
at a solution point, the Mach number, static temperature, velocity magnitude, 
stagnation pressure, and stagnation temperature are also printed. The mass 
flow rate across every solution plane, calculated by trapezoidal rule inte- 
gration in subroutine MASS, is also printed. 

The initial data and subsequent solution surfaces for the boundary layer 
calculation are printed by subroutine DEC0DE. The. output consists of position, 
velocity components, pressure, density, temperature, velocity magnitude, stag- 
nation pressure, stagnation temperature, bleed rate, and certain integral para- 
meters. 
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SECTION III 

SUBROUTINE DESCRIPTIONS 


1 . INTRODUCTION 

In this section, a brief description is given of the function of each sub-^ 
routine in the computer program. This information supplements the information 
available in the form of comment statements within the program. 


2. 0VERLAY (0,0) 

LINKOO . This program routine is the main control routine in 0VERLAY (0,0), 
the resident overlay. LINKOO first calls 0VERLAY (1,0) for data input, para- 
meter initialization, and, if desired, internal generation of the initial data. 
LINKOO then calls 0VERLAY (3,0) and 0VERLAY (4,0) to perform the flow field in- 
tegration. Most of the program constants and input parameters have their default 
values specified in LINKOO. Moreover, the reader and printer call numbers, de- 
noted by IRE and IWR, respectively, are initialized in LINKOO. 

B0UNDY . This subroutine is used for the specification of both the fore- 
body/centerbody and cowl geometries. The version of B0UNDY supplied with the 
program assumes that both the forebody/centerbody and the cowl are axi symmetric. 
More arbitrary geometries can be described by suitably modifying or replacing 
this subroutine. If B0UNDY is replaced, the subroutine argument list must be 
Identical to the existing list. The parameters in the argument list are de- 
fined as follows. 

XABS x-coordinate 

Y y-coordinate 

Z z-coordinate 

ALPHA polar angle defined by tan“^ (z/y) 

RB0DY radius of either forebody/centerbody or cowl 

BNX x-component of outward body normal unit vector to forebody/ 
centerbody or cowl 

BNY y-component of outward body normal unit vector to forebody/ 
centerbody or cowl 

BNZ z-component of outward body norma) unit vector to forebody/ 
centerbody or cow.1 

J If J=l, forebody/centerbody geometry is to be specified. 

If J=2, cowl geometry is to be specified. 
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K 


If K*0, read In XABS and ALPHA, and compute RB0DY and the cor- 
responding Y.and I. If K**!, read in XABS and ALPHA, and compute 
RB0DY, Y, Z, BNX, BNY, and BNZ. If K*2, read in XABS., Y, and Z 
(coordinates of a point not on the. body), and compute RB0DY, 
ALPHA, Y, Z, BNX, BNY, and BNZ at the point on the body where 
the projection of the body normal in the (y,zj-piane passes 
through the originally specified point. 


IH. E R.MP; This subroutine computes the temperature, sonic speed, and the 
thermodynamic parameter The assumed thermodynamic model is that of a 
thermally and calor.ically perfect gas. Other thermodynamic models may be in- 
corporated into the analysis by suitably modifying orreplacing this subroutine. 
If subroutine THERM0 is replaced, .the argument list must be identical to the 
existing list. The parameters in the argument list are defined as follows. 


P pressure (p) 

R0 density (p) 

T" temperature (t) 

A sonic sp.eed (a). 

C0EFF 5=(l/pt)(ap/as)p, where s is the entropy per unit mass 

KT If KT=1, compute. T. If KT?«1 , do not compute , T. 

KA If KA»1 , compute A'. If KA/1 , do not compute A. 

KC If KC=1 , compute C0EFF. If KCj^l , do not compute C0EFF. 

Thermodynamic property information is supplied to subroutine THERMj!) through the 
named common block GASl. This common block contains the following parameters: 
RUNiy (universal qas constant), TDT (temperature data array), WTM0L (molecular 
weight data array), CP (specific heat data array), NGASl .(an Integer to denote 
the number of elements in. the TDT, WTM0L, and CP arrays), and KGASl (an Integer 
to denote which gas model Is to be used). These parameters are Included for 
modifying subroutine THERM0 to, include real gas effects. 


stag . This subroutine computes the Mach number, the stagnation pressure, 
and the stagnation temperature. The assumed thermodynamic model is that of a 
thermally and calorlcally perfect gas. Other thermodynamic models may be in- 
corporated into the analysis by suitably modifying or replacing this subroutine. 
If STAG is replaced, the argument list must be identical to the existing list. 
The parameter.s in the argument list are defined as follows. 


P static pressure 

R0 . static density. 

T static temperature 

A son.ic .speed. . 

Q velocity magnitude 

XM — Mach number 


PT 

TT 


stagnation pressure- 
stagnation temperature 


STAG through the 

specIfTc^thafpy! 'Any°two"vartaburm^''^h'’''®“‘'''®^‘'®''®*‘^’ ‘<h"P«>'ature, and 
The assumed thermodynamic model is computed given the remaining two. 
gas. other thermSmic mSde s calor.cally perfect 
ably modifying or replacInrthi^suLouti^e analysis by sult- 
the argument list niust be fdent cai Is replaced, 
the argument list are defined as follows. ^”^"9 parameters In 


P 

R0 

T 

H 

K0PT 


pressure 
density 
temperature 

specific static enthalpy 

If !JSdt°J* R0 and T given P and H. 

If T and H given P and R0. 

If ^Sdt’I* R0 and H given P and T. 

If ^ given P and T. 

If K0PT-5, compute R0 and T given H and P. 

SlTpS’cloThlJJj: JXsr^‘'°" tuhroutlnoJHERHH through 

using an analytt^^l^'^diff^rent1Sd°forIn^f''thI°fhl'’’^’“‘® temperature by 
that temperature derivatives can be exor^ccL^f^ thermal equation of state so 
pressure and density. ThrLsSlS theEm^^ derivatives of 

perfect gas. Other thermodynamic models that of a thermally 

by suitably modifying or reolacinn thi'c cuk'^ bs incorporated into the analysis 

argument list mist SI iSentfJal SI Shi ellsSSlS fS;t 

argument list are defined as follows. ^ ^ ^^® Parameters in the 


P 

R0 

DPY 

OR0Y 

DTY 


pressure 
density 

first partial -spatial derivative of pressure 
first partial spatial derivative of density 
first partial spatial .derivative of temperature 


nteTSI'bfS?S'Sil!"'°™“'“" through the 
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dient^ho 9rd- 

Ldw ( 6 ). othei* transport propmJ forinu?at(ons S'*®" '=d' Sutherland' 

ysis by suitab-ly modifyin/or rlolaci^rthu J k ^nal' 

the argument list must be subroutine, if lvIS Is replaced. 

argument list are defined as follows. existing list. The parameters in the 


T 

P 

VIS 
DTY 
DPY 
DVISY 
K0PT 


static temperature 
static pressure 
dynamic viscosity 

first partial spatial derivative of temperature 
first partial spatial derivative of pressure 
first partial spatial derivative of. viscosity 
If K0PT=1 , compute VIS. 

If K0PT=2, compute, DVISY. 


common block^GASzf ^The*"Smon block^GAS2^contai^'"°jr'^f^ through the named 
TDL (temperature data arraT VISD (v?sLsU^h follow ng parameters: 

tivity data array), NGAS2 (an intpaprtn°wI^'^^.^®J? a>"»"ay), C0ND (thermal conduc- 
TDL, VISD, and C0ND arrays) and^KGA^?^?;,n the number of elements in the 
property model iftn hi -integer to denote which transport 

LVIS tTusHabulaJ parameters are included for modifying 

quadrf^^•n^e^po?ati^^^^ the Ihree coefficients in the univariate 

equations. The^system^is'^solvpd^^*?i^°'"r^°^'^^-'^^ ^ System of. simultaneous linear 
pivoting. This s^b^out L i ?rm ? "^^‘"’i^ation with complete 

with ri^etric^ioefffcle^ \hrsvstPm linear equations 

el1n.1natl0h with pivot „g in Jh^min “=*"3 Gaussian 

version of: IBM library Lb™utine. Gas^ ®^^^ .. a mdlfled 

meters fo*r 'use^1n*the*’Sp.utat?on*'™' the, wall transpiration/suction para- 

the bWl’ven^the .ltSlse cu^lllneS? * Point on 

3. -0VERLAY_(1.O) 

All 1ff^ara™Mrs™?d'”nUl^'"dauire entered'’lh™LINK!lo''’*l^ ‘’-j’' 

enter this Information. After the iJIpu? dat^have bee7en‘ttel!'sJbro“5t1ne° 
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INITIL is called 'for testing for input errors and performing parameter initial* 
ization. If the flow property field on the initial-value plane is to be inter- 
nally generated, the appropriate subroutines are then called from subroutine 
INITIL. After subroutine INITIL has been called, preliminary output is. printed 
by LINKIO. 

INITIL . This subroutine is called from routine LINKIO, and is. used to test 
for errors on selected input parameters and to perform paramete.r initialization. 
If an input parameter error is found, an appropriate message is generated and 
the program execution is aborted. After the selected input parameters have 
been tested, parameter initialization is performed. If the flow property field 
on the initial-value plane is to be internally generated, the apprppriate sub- 
routines are then called. 


4. PVERLAY (1,1) 

LINK!! . This program routine is the main control routine in OVERLAY (1,1), 
and is used in conjunction with subroutines TMC0NE, INT6RT, L0GIC, C0NVR6, 

0DESLV, and FUNCTM to internally generate the supersonic flow property field 
on the initial-value plane if the forebody is conical ahead of the initial- 
value plane. The internally generated initial-value plane is obtained by an 
approximate technique which employs the Taylor-Maccoll solution for the flow 
about a circular cone at zero incidence. A superposition procedure is used to 
obtain. an approximation to the flow about a circular cone at nonzero angle of 
attack by neglecting the cross flow effects. This superposition procedure 
effectively amounts to computing the flow turniag angle in the meridional plane 
of the given solution point, and then obtaining the flnw properties at that 
point by applying the Taylor-Maccoll solution for a cone half-angle equal to 
the flow turning angle. The shock wave angle is then measured from the origi- 
nal streamline direction in the appropriate meridional plane- It must be em- 
phasized that this is only an approximate technique, giving the well accepted 
Taylor-Maccoll solution at zero angle of attack, but becoming increasingly less 
accurate as the angle of attack is increased. 

TMC0NE . This subroutine is one of the subroutines used for the internal 
generation of the initial-value plane. TMC0NE is called from routine LINKll 
and performs initialization and post calculation storage of the Taylor-Maccoll 
solution used in computing the initial-value plane supersonic flow property 
field. Ths assumed thermodynamic model is that of a thermally and calorically 
perfect gas. Other thermodynamic models may be incorporated into the analysis 
by suitably modifying this subroutine- Thermodynamic property information is 
supplied to subroutine TMC0NE' through the named, common block OASl. 

INTGRT. This subroutine is one of the subroutines used in internally gener- 
ating the supersonic flow property field on the initial-value plane. INTGRT in- 
tegrates the Taylor-Maccoll equation inward towards the. forebody until a point 
is reached where the normal component of velocity is sufficiently small. 


18 


,LP6IC.. This subroutine is one of the subroutines used in internally gen- 
erating the supersonic flow property field on the initial-value palne. Lj5GIC 
is used to control the computation of the Taylor-Maccoll solution in a particu- 
lar meridional plane. The assumed thermodynamic model is that of a theinnally 
and calorically perfect gas. Other thermodynamic ’models may be incorporated 
into the analysis by suitably modifying this subroutine. Thermodynamic property 
information is supplied to subroutine LJ5GIC‘ through the named common block GASl . 

CgNVRG. This subroutine is one of the subroutines used in internally gen-, 
erating the supersonic flow property field on the initial-value plane. CPNVRG 
■js.used in testing if the normal component of velocity has vanished at the body 
in the external flow field integration. . 

£D|SL\/. This subroutine is one of the subroutines used in internally 
generating the supersonic flow property field on the initial-value plane. 0DESLV 
integrates the Taylor-MaccolT equa^on by employing a fourth-order Runge-Kutta 
method. ’ 


This function is used in internally generating the supersonic flow 
property field on the initial-value plane. FUNCTM evaluates the Taylor-Maccoll 
equation. The assumed thermodynamic model is that of a thermally and calorically 
perfect gas. Other thermodynamic models may be incorporated into the analysis 

this function. Thermodynamic property information is sup- 
plied to FUNCTM through the named common block GASl. 


5. 0VERLAY (1,2) 

This program routine is the main control routine in 0VERLAY (1,2) 
and IS used in conjunction with subroutines STENCL, REFLK, and BI in converting 
the output of the Jones algorithm into a format which is compatible with the 
bicharacteristic algorithm. The Jones algorithm is contained in 0VERLAY (2,0) 
and IS used to internally generate the supersonic flow property field on the 
initial-value plane. The assumed thermodynamic model employed in LINK12 is that 
of a thermally and calorically perfect gas. Other thermodynamic models may be 
Incorporated into the analysis by suitably modifying this routine. Thermodynamic 
property information is supplied to routine LINK! 2 through the named common block . 
GAdi . 


SIENCL. This subroutine is one of the subroutines used for conversion of 
the Jones algorithm output so that it can be used as initial data for the super- 
sonic flow solution. STENCL determines the stencil of field points to be used 
in the bivariate interpolation of the initial data. 

REFLK. This subroutine is one of the subroutines used for conversion of 
the Jones algorithm output so that it can be used as initial data for the super- 
sohic flow solution. _ REELK performs field point reflections about the plane 
of flowLaymmetry. 
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BI. This subroutine is one of the subroutines used for- conversion, of the 
Jones algorithm output so that it. can be used as initial data for the Supersonic 
flow solution. BI determines Jthe coefficients. Of the least squares quadratic 
bivariate polynomial used for interpolation of the initial data. . The system 
normal equations which determines the polynomial coefficients is solved by calling 
subroutine SYMSIM. . 


6. 0V£RLAY (1,3) 

-LINK13 . This program routine is the main control routine in (JVERLAY (1,3) 
and is used in conjunction with subroutines SADAMS, PEQS0, INTER, C(8EFF, BINQ, 
and VEL to internally generate the initial data for the forebody/centerbody boun- 
dary layer computation. The initial data generation is accomplished by using 
the Adams finite-difference algorithm (5), which is incorporated, into the compu- 
ter program .in OVERLAY (1,3). The Adams algorithm generates the initial data at 
a specified axial, location, assuming th't the forebody/centerbody geometry is 
conical ahead of that location. It is also assumed in the analysis that the 
wall temperature is a specified constant, that the thermodynamic model. Is that 
of.a thermally and calorically perfect gas, and that dynamic viscosity is given 
by Sutherland's Law. (6), Other thermodynamic and molecualr transport models 
may be incorporated into the analysis by suitably modifying the subroutines in 
OVERLAY (1,3). Thermodynamic and molecular transport property information is 
supplied- to routine LINKl 3 through the named common blocks BASl and GAS2. 

SADAMS. This subroutine is one of the subroutines used for the internal 
generation of the forebody/centerbody boundary layer initial data. SADAMS is 
the main control routine in the Adams finite-difference algorithm. 

PEQS0 . This subroutine is one of the subroutines used for the internal 
generation of the forebody/centerbody boundary layer initial data. PEQS0 solves 
the system of parabolic differential equations generated in the analysis. 

inter . This, subroutine is one of the subroutines used for the internal 
generation of the forebody/centerbody boundary layer initial data. INTER '■> 
used to .perform Lagrangi an interpolation. 

C0EFF . This subroutine is one of the subroutines used for the internal 
generation of the forebody/centerbody boundary layer initial data. C0EFF is 
used to obtain the .surface properties given the surface pressure distribution. 

BIND . This, subroutine is one of. the subroutines used for -the internal 
generation of the forebody/centerbody boundary layer initial data. BINQ solves 
the system of simultaneous linear equations used in the Newton iteration scheme. 

VEL . This subroutine is. one of the subroutines used for the internal gen- 
eration of the forebody/centerbody boundary layer initial data. VEL is used to 
interpolate the surface properties determined by subroutine C0EFF. 
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7. OVERLAY (1,4) 


LINK14. This program routine is the main control routine in 0VERLAY^O«4) 
and is used to initialize the transformed variables used in the boundary layer 
computations. 


8. OVERLAY (2,0) 

LINK20 . This program routine is the main control routine in 0VERLAY (2,0) 
and is used in conjunction with subroutines JJ8NALG, FUNClli, CHI, FA, JSHPCK, 
SRFACE, PMC ZER0SN, FAl , FN, RKINIT, RKSTEP, DWEG, MAINM, MINIV, INVERT, DS0LV, 
and B0DY to internally generate the supersonic property field on the initial- 
value plane if the forebody is conical ahead of the initial-value plane. The 
internally generated initial-value plane is obtained by using the Jones finite 
difference algorithm (4). The Jones algorithm has been incorporated into the 
computer program in 0VERLAY (2,0). It is assumed that the thermodynamic model 
is that of a thermally and calorically perfect gas. Other thermodynamic models 
may be incorporated into the analysis by suitably modifying the subroutines in 
0VERLAY (2,0). Thermodynamic property information is supplied to routine LINK^JU 
through the. named common block GASl. All of the subroutines in WLAY (2,0) . 
are used for the internal generation of the supersonic flow initial data. 


JgNALG . J0NALG is the. main control routine for. the Jones algorithm. 

FUNCTN. FUNCTN computes the surface residual mass flow functions which 
are to be minimized. 


CHI . CHI computes certain output parameters. 

FA . FA computes gradient parameters that are used in subroutine PMC. 

jSHgCK . JSH0CK computes the flow properties immediately downstream of the 
bow shock wave. 

SRFACE . SRFACE finds the initial approximation to the surface solution. 

pmC . PMC is used for numerically integrating the governing differential 
equations .. 

ZERgSN* ZER0SN is_used to calculate the solution for a circular cone at 
zero incidence. 

FAl. FAl computes the derivatives used in the Runge-Kutta integration for 
the zero incidence solution. 

FN.. FN controls the Runge-Kutta integration for the Zero incidence solution. 

RKINIT. RKINIT Computes functions used in the zero incidence solution. 


RKSTEP . RKSTEP -performs one Run^e-Kutta integration step for the 2ero in- 
cidrtece solution. 

DWE6 . _PWE6 is used for computing zeros of an equation. 

MAINM . MAINM performs functional minimization applied for obtaining a 
minimum normal velocity at the body surface. 

MINIV. MINIV is used for functional minimization. 

INVERT . INVERT is used to perform matrix inversion. 

DS01V . DS0LV uses Gaussian elimination for solving systems of linear 
equations. 

SgDY . B0DY is used to define conical geometries. 


9. OVERLAY (3,0) 

LINK30 . This program routine is the main control routine in OVERLAY (3,0). 
LINK30 calls the supersonic flow integration control overlays, 0VERLAY (3,1), 
0VERLAY (3,2), and ^yVERLAY (3,3). If a program restart is specified, the last 
secondary level overlay that was in use at the time the initial execution was 
terminated is called. 

REDIST . This subroutine is used to perform point redistribution on the 
solution plane at the axial station of the cowl lip. This point redistribution 
is required to obtain a uniform point distribution and to obtain streamlines 
which lie in the stream surface formed by the cowl boundary. The redistributed 
points are arranged symmetrically in the computed sector. These points lie on 
rays which have equal angular increments from one another, with the points on 
each ray being spaced at equal radial increments. The flow properties at these 
points are obtained by bivariate interpolation. If the viscous flow option is 
specified and the internal flow field integration option in which shock waves 
are not discretely fitted is specified, point redistribution is also performed 
on the forebody flow field solution plane immediately upstream of the solution 
plane. located at the cowl lip axial station. 

STNSEL . This subroutine is employed to determine the solution point on a 
given .solution plane that is nearest to an arbitrary point whose coordinates 
are supplied in the argument list. After the proper solution point has been 
determined by a search of the computed sector, it is used as the base point 
of the nine point stencil used in formulating the quadratic bivariate inter- 
polation polynomial. STNSEL calls subroutine BIPTS in the course of determin- 
ing the proper solution point and in formulating the interpolation polynomial 
coefficients. 

REFLEK . This subroutine performs reflections, about the Coordinate axes, 
of the solution points in the computed sector on a given solution plane when 
flow symmetry exists. The solution point reflections are required to obtain 
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the appropriate fit point stencils used in forinulating the univariate, bi- 
variate, and trivariate interpolation polynomials. REFLEK is used to reflect 
points for tho case of one -plane of flow symmetry, the case of two planes of 
flow symmetry, and the axi symmetric- flow case. Both streamline and shock wave 
points are reflected in REFLEK. 

.! XSTEP. This subroutine determines, the axial-marching. step allowed by the 
Courant^Friedrichs-Lewy (CFL) stability criterion. Except in the vicinity of 
an internal -Shock. wave intersection with a solid boundary, the marching step 
computed by XSTEP is that used to regulate the distance between successive 
solution planes. In the vicinity of an internal shock wave-solid boundary inter- 
section, special constraints are used to determine the marching step. XSTEP 
applies the CFL Stability criterion to all streamline points in the computed 
sector. The acutal marching step is taken as the allowable x-step at the most 
restrictive point. The stability criterion is not applied to the shock wave 
points. Moreover, the internal flow field shock wave points are ignored in de- 
fining the convex hull of the finite difference network when application of the 
stability criterion is made to a streamline point. 

MASS . This subroutine is used to compute the mass flow rate across a given 
solution plane. The mass flow rate is numerically ascertained by a trapetoidal 
rule integration procedure in which the incremental mass flow rate is computed 
through an elemental triangle formed by three adjacent points. The sum of these 
incremental mass flow rates is the total mass flow rate. If flow symmetry exists, 
on'y the mass flow rate in the computed sector is found by integration. The mass 
flow rate across the entire plane is then obtained by use of an appropriate multi- 
plier. Special logic is used to trace the internal shock-wave in the integration 
procedure when the mass flow rate across an internal flow field solution plane is 
being computed. 

PRN0UT. This subroutine is used to print the supersonic flow integration^ 
results for all solution planes. Both the external flow field and the internal 
flow field solution planes are printed by PRN0UT. In addition, the initial- 
value plane, the redistributed data plane at the cowl lip axial station, and 
the restart plane are output by PRN0UT. Three print options (specified in the 
input by the input, parameter KPRINT) are available: body streamline points and 

shock wave points, all solution points, and all solution points and shock wave 
parameters (incident normal Mach number and shock wave surface normal unit vector 
components). The solution.points along the space curves formed by tb.'». intersections 
of the internal shock wave with the solid boundari.es are not printed by subroutine 
PRN0jJJ. They are printed by Subroutine .SHKRFL. 

UNlpTr This subroutine selects the fit points used in formulating the 
.in I vaFi ate interpolation polynomials that are used to describe.the shock wave 
radius and the shock wave angle along the Curve defined by the intersection of 
the shock wave, with agiven solution plane. Three adjacent shock wave solution 
points constitute the fit point stencil in regions that are away from an internal 
shock wave-solid boundary intersection. In the region of an intersection, the 
fit point stencil may be appropriately expanded to be in accord with the Courant- 
FriedriChs-Lewy (CFL) stability criterion. 
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UNI FIT . This subroutine is used to determine 

rati c"unTOri ate which are supplied to UNIFIT 

nomial are determined by a to tjree data pom 

Srk-sllU“ th.e_ coefficients. 

BIPTS . This subroutine selects the fit. poin^ forTintraplanar flow 

cuadratTTbivariate ^oterpal at ion polynomials that are se ^bors const!- 

property determination. A base point and its eight imm^ ^ 

sL':L'1s.k s:t “-F> 

Es^riKS.'s "!.« »"s;,r.s: rsan...!.' »• 

Shockwave. . *u 

BIFIT. This subroutine is ^^^®^®Je^stx^coeff^ polyrto- 

ratic'bf^riate interpolation polynomial. The s x o t supplied 

m?ri are determined by a least squares fit of nine^dat^ 

to BIFIT through the earned common^ b^^^ Rations (normal equations) which de- 
fe'^ln'ertSf *c^f“4cnrrnT4 This syste™ of 1 ,n 

equations has a symmetric coefficient matrix. 

TRIPTS. This subroutine selects the fit points used i^n^fo^uUt^ 
the lineaTtri variate interpolation Po]yj°"’'’®^ fSmulating the linear tri- 

polation polynomial. solution points are used in formulating 

variate interpolation polynomial. Fourteen h polynomial is used 

Th^quadratiC trivariate interpolation shock Jiave surface 

rt?ersrdrS/?’hrsK^'K r?hs fo™a. a S 0 l,h 

boundary. 


linear 


TRIFTU. This subroutine is used to determine the polynomial 
trivarTatT'interpolation polynomial, /^he four c to TRIFTU through 

Trl delermined by a fit. to four Oe^Erc' s'raireS'ti solve the system of 

TRIFTD. This subroutine !{y®Q ^°l^®^Xhe'Oight^ i^ this poly- 

supplIPl to TRino through .the mHlins formal equatlohs) 

£iJ^a*e?e^rSlres‘5?e^?rerfllIe5;?“’H"hra“UrerSfi;^ has a sy-etr.c co- 
efficient matrix. .. . 

SHKINT. This, subroutine is used to evaluate^the JJJ^g^OadiuS and the' shock 
polynSHad that are by the Intersection of the shock, wave 

k rfltej’ro?u‘t^on'“pra^e!''1h:"l2irrr^ polynomial coefficients are 
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supplied to SHKINT through the named common block SINTRP. The Independent vari- 
able Is the pOUr. angle THETA, which Is transmitted through the subroutine call 
statement . 

PRPINT. This subroutine is used to evaluate the bivariate and tri variate 
interpolation polynomials that are used for flow property determination on solu- 
tion planes, shock wav£s, and solid boundary stream surfaces. The interpolation 
polynomial coefficients are supplied to PRPINT through the named common block 
INTRP. The independent variables are the three position coordinates XA, Y, and 
Z, which are transmitted through the subroutine call statement. 

FDERIV . This subroutine is used to compute the first partial derivatives 
of the flow properties on a solution plane, a shock wave, or a solid boundary 
stream surface. The derivatives are obtained by evaluating an analytically dif- 
ferentiated form of the appropriate interpolation, polynomial . The Interpolation 
polynomial coefficients are transmitted to FDERIV throagh the named common block 
INTRP. The independent variables are the three position coordinates XA,. Y, and 
Z, which are transmitted through the subroutine call statement. 

VECT0R . This subroutine is used in conjunction with subroutine SH0CK to 
determine the components of the unit vector ^ which is used in the parameteriza- 
tion of the wave surface compatibility relation. The unit vector *0 is orthogonal 
to the velocity vector that is downstream of the. shock wave . at the shock wave so- 
lution point, and has its projection on the (y,z)-plane lie in the meridional 
plane which passes through the shock wave solution point. 

PLANAR . This subroutine is used to compute the parameters employed in the 
formulation that represents the shock wave surface. This formulation is then 
Used in determining the intersection point of either a streamline or a bichar- ■ 
acteristic with the shock wave surface. Additionally, PLANAR initializes some 
parameters used in determining the intersection point of a bi characteristic with 
the stream surface formed by a solid boundary. 

INTSCT. This subroutine is used to compute the intersection point of either 
a streamline with the shock wave surface, or the intersection point of a bi char- 
acteristic: with either the shock wave surface or the stream surface formed by a 
solid boundary. The determination of the intersection point coordinate? 1: per- 
formed In an iteration loop which uses the secant method to relax the difference 
in the radius of the point of intersection obtained from integration.o.f the equa- 
tion for a streamline or bi characteristic and that obtained from the -appropriate 
surface formulation. 

FUNCT. This subroutine is used in conjunction- with subroutine INTSCT to 
determine a streamline-shock wave Intersection point, or-the intersection point 
of a bi characteristic with either a shock wave or a solid boundary. The axial 
position of 'the assumed intersection point is supplied to FUNCT through the 
subroutine call -statement. From the given axial position, FUNCT computes the 
corresponding y and z Coordinates of the assumed intersection point by both 
integrating the equation for a streamline or bicharacteristic, and by evaluating 
the appropriate surface formulation. The difference in the radius obtained from 
the streamline or biCharaCteristic equation integration and that obtained from 
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the surface formulation Is reduced to within a specified tolerance of zero by 
the Iteration method used In subrautine INtSCT^. 


RADIUS , This subroutine Is used to compute the radius of a point .on the 
shock“wave surface or on the stream surface formed by a solid boundary. The 
axial position and the polar angle of the point are supplied to RADIUS through 
the subroutine call statement. To obtain the shock wave radius at the desired 
point, a linear Interpolation is performed in the meridional plane of the point 
between two space-curves which are defined by either a shock-wave solution plane 
Intersection or by a shock wave-solid boundary Intersection. To obtain the body 
radius at the desired point, subroutine B0UNDY is called, 

XFIT . . This subroutine Is used to curve fit, as a function of the polar 
angle, the axial (x) position of an Internal shock wave-solid boundary Inter- 
section. A quadratic polynomial expressed.ln terms of the polar angle is used 
for this representation, 

FRCFNS . This subroutine Is -used to compute the molecular transport forcing 
terms used In both the governing equations end the compatibility relations. 
FRCFNS-is called by subroutines S0LVE and SHjSCK. The molecular, transport terms 
can be Included In the computation of the external flow field about the forebody 
or In the computation of the Internal flow field In which shock waves are not 
discretely fitted. The program option In which Internal shock waves are dis- 
cretely fitted does not have the capability to Include the Influence of molecu- 
lar-diffusion In the computation, but rather assumes the flow to be Inviscid 
and adiabatic. 

SgLN . This subroutine calls subroutines L0CATE and S0LVE for computing 
either a solid Jjoundary solution point or an interior solution point. 

L0CATE . This subroutine is used to compute the locations of and the flow 
properties at the streamline and bicharacteristic Intersection points for the 
standard Interior point, standard solid boundary point, shock-modified Interior 
point, and shock-modified solid boundary point unit processes. The point loca- 
tions and flow properties at these points are transmitted to subroutine S0LVE 
through the named common blocks RELAYl and RELAY2. Subroutine S0LVE then solves 
the system of compatibility relations to obtain the position of and the flow 
properties at the solution point. 

^ S0LVE . This subroutine Is used to solve the system of compatibility eqya-_ 

I tions for the standard Interior point, standard solid boundary point, shock- 

: modified Interior point, and shock-modified solid boundary point unit processes. 

If the viscous and thermal diffusion terms are to be Included In the computation, 

I S0LVE calls subroutine FRCFNS. The system of five compatibility relations Is 

solved by. calling subroutine GELG.- 

I SH0CK . This subroutiJie Is used to compute the solution for all field-shock 

F wave points. For the bow shock wave points, the free-stream flow conditions are 

I used for the upstream flow properties. For.the Internal flow field shock wave 

! points, an Interior or sol Id boundary point unit process Is applied to obtain 

I the upstream flow properties. Body solution points on the downstream side of 
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the cowl -lip shock wave or oh the downstream side of a reflected Internal shock 
wave are computed using the solid boundary~shOck wave point unit process (sub- 
routine BSH0CK). 

JUMP. This subroutine employs the Rankine-Hugoniot relations to compute 
the downstream flow properties at a shock wave solution point. The shock wave 
surface normal unit vector components and the flow properties upstream of the 
shock wave are entered in the argument list. The computed downstream flow 
properties are also transmitted in the argument list. The assumed thermodynamic 
model is that of a thermally and calorically perfect gas. Other thermodynamic 
models may be incorporated into the analysis by suitably modifying this sub- 
routine. Thermodynamic property information is supplied to subroutine JUMP 
through the named common block GASl. 

TGRAD. This subroutine computes the spatial gradients of temperature by 
using an analytically differentiated form of the thermal equation of state so 
that temperature derivatives can be expressed in terms of the derivatives of 
pressure and density. The assumed thermodynamic model is that of a thermally 
perfect gas. Other thermodynamic models may be incorporated into the analysis 
by suitably modifying or replacing this subroutine. If TGRAD is replaced, the 
argument list must be identical to the existing list. The parameters in the 
argument list are defined as. follows. 


T 

P 

R0 

DPX.DPY'yDPZ 

DPXX,DPYY,DPZZ 

DR0X,DRj5Y,DR0Z 

DR0XX,DR0YY,DRPZZ 

DTX,DTY,DTZ 

DTXX,DTYY,DTZZ 


temperature 

pressure 

density 

first partial derivatives of pressure with respect 
to X, y, and z, respectively 

second partial derivatives of pressure with respect 
to X, y, and z, respectively 

first partial derivatives of_density with respect 
to X, y, and z, respectively 

second partial derivatives of density with respect 
to.Xj.y, and z, respectively 

first partial- derivatives of temperature with 
respect to x, y, and z, respectively 

second partial derivatives of temperature with 
respect to x, y, and z, respectively 


Thermodynamic property information is supplied to subroutine TGRAD through the . 
the named common block GASl. 


TPR0PS. This Subroutine computes the dynamic viscosity, 
ductivity, and their Spatial gradieats. The assumed functional form for dynamic 
viscosity is given by Sutherland's Law (6). The assumed functional form for 
thermal conductivity is expressed in terms of the dynamic viscosity and laminar 
Prandtl number. Other transport property formulations may be incorporated into 
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suitably modifying or replacing this -subroutine. If TPR0PS is 
replaced, the argument list must be Identical to thp eyicf<nn nef ~ 'S 
meters In the argument list are defined as follows ® 


T 

P 

VIS 

C0N 

DTX,DTY,DTZ 
DPX.DPY.DPZ , 
DVISX.DVISY.DVISZ 
DC0NX,DC0NY,DC0NZ 


Static temperature 
static pressure 
dynamic viscosity 
thermal conductivity 

first partial derivatives of temperature with 
respect to x, y, and z, resp.ectively 

first partial derivatives of pressure with respect 
to X, y, and z, respectively 


first partial derivatives of viscosity with respect 
to X, y, and z, respe.ctively ^ 


first partial derivatives of thermal conductivity 
with respect to x, y, and z, respectively 


Transport property Information 1 
named common block GAS2. 


s supplied 


to subroutine 


TPR0PS through the 


10. 0VERLAY (3.1) 

program routine contains the control logic used in the conmu 
tation of the supersonic external flow field about the forebodv LINK31 ic 

MKi ■■ ra..s"s.r.;“ 

current solution plane and th^l^sE^Slutlon'planrwher^^ SE®de 

Doints If outermost ring of existingjnterior field 

U successive point additions, a specified number 

reached, point deletion Is performed. Here, selected 

i«rIointr „d'? 2‘Lr amy, wh™e the bo™ sSIcf 

wave points and the. body streamline points are retained. 


11. 0VERLAY (3,2)_ 

. seconde^^"?e1el%5S5,;;’l5?S‘f?’c;?^:J"?,:: f^ERLAY (3.3J._ ,VERLAY.(3.2) 
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BSH0CK. This subroutine is used to compute^ the flow 
on the body that is downstream of either the cowl lip shock wave or an internal 
reflected shock wave (solid body-shock wave point unit process). This, sub- 
routine is used only in the option where the internal flow field integration in 
which shock waves are discretely fitted is employed. 

SPRJCT. This subroutine is used in the internal flow 
project the Internal shock wave from the current initial-value, plane to the cur- 
rent solution plane. The projected intersection of the internal shock wave with 
the solution plane is then used to determine which streamlines do and do not 
penetrate the shock wave. 

PENTRE. This subroutine is used in the internal flow field integration to 
control the computation of an interior field point when the streamline has P®|^®" 
trated the internal shock wave. If the streamline-shock wave intersection point 
is sufficiently close to the current solution plane, an interior point unit pro- 
cess on the downstream side of the shock wave is not performed. Instead, in 
this case, a streamline projection onto the solution plane and subsequent 
property interpolation in this plane is performed. /^t®»'^®^Jvely, if the strea^^ 
line-shock wave intersection point is sufficiently far from ^ 

an interior point unit process is performed on the downstream side of the shock 

wave. 

TRACE. This subroutine is used in conjunction with subroutine INTSCT 
compute the intersection point of a streamline with an internal shock wave. With 
the intersection point coordinates detemined, the trivariate interpolation poly- 
nomials are evaluated to obtain the flow properties at the intersection point. 

S0LINT. This subroutine is used in the internal flow field integration to 
determine the intersection point coordinates of and the flow properties at an 
interior field point when the streamline penetrates an internal shock wave with 
the intersection point being sufficiently close to the current solution plane, 
S0LINT is called from subroutine PENTRE. 

SWITCH. This subroutine is used in the internal flow field integration to 
perform the post computation interchange of indices when a ^^rnn«p'" 

tially appeared to Intersect the internal shock wave ultimately 
quently, a standard interior point unit process is used for this J."stead 

of a shock-modified interior point unit process performed on the downstream 
side of the shock wave. 

SHKRFL. This subroutine contains the control logic used for calculating 
an internal shock wave-solid boundary intersection. 

INTRFL. This subroutine performs the necessary point 
symmelryi^istrfor the solutiSn points oa the space curve defined by the inter- 
section of the internal- shock wave with a solid boundary. 

ISH0CK. This subroutine is used in the internal flow 
compute the incident shock wave upstream and downstream flow properties at a 
point where the incident shock wave intersects a solid boundary. 
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ProDeflfe^at thl ^ compute the position of and the flow 

defined bv ^ streamline with the space curve 

fined by ihe intersection of the incident shock wave with a_solid boundary. 

This subroutlns Is us6d with subroutinA +a AAmAti 4 .r^ 4 .ux 

flow properties on the body downstree/rf er?n™e™i ?e'??ecte3 

12. 0VERLAY (3,3) 

tatlo^e IIlJlrsl1c“nt™rn‘If%fo“w‘?leld‘?r 

(3.3) is a secondary level overlay which is called from routine LINO^ 

13. OVERLAY (4,0) 

LINK4^S^irJhe'^boSnXrUye^ J^pStatiln o^bD?!! OVERLAY (4,0). 

and the cowl LiNKan raiic fl on both the forebody/centerbody 

U? (S^rWRuf'Ur.'Lf^VERwlJ^ar^ ^VER- 

ture “offfcleots and geodesic curva- 

a^fdefined al filfows^ P^^^n^^ters in the argument list 


KB0DY 

KPLANE 

XBASE I 
XBASES 


If KB0DY=1, forebody/centerbody geometry is to be 

S ps Cl T 1 6 Q • 

If KB0DY=2, cowl geometry is to be supplied. 

If KPLANE=1 , axial location is for current initial- 
value plane. 

axial location is for current solution 

plan 6* 

Base axial coordinate of the .current initial-value 
plane. 

Base axial coordinate of the current solution plane. 


base Si„.2!' Thl™e)il?p7f'’xcSRv‘sSo!l?l^*’’4Sn[ 
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KBjKDY If KB0DYal , forebody/centerbody_ geometry is to be 

specified. 

If KB0DY=2, cowl geometry is. to be specified. 

XBASEI Base axial coordiante of the current initial~value 

plane. 

XBASES Base axial coordiante of the Current solution plane. 

EINTRP . This subroutine determines the external flow properties for the 
boundary Tayer computation by interpolation of the supersonic core flow solution. 

WALL . This subroutine calculates the wall condition parameters for compu- 
tation of the boundary layer. 

WTEMP . This subroutine determines the wall temperature boundary condition 
for use in the boundary layer computation. Both temperature and normal tempera-: 
ture gradient boundary conditions can be specified. Subroutine WTEMP is called 
from subroutine WALL. 

TURB . This subroutine computes the turbulent eddy viscosity and turbulent 
eddy thk-.''.al conductivity for use in the boundary layer computation. The version 
of TURB supplied with the program employs an isotropic two-layer eddy diffusitivity 
model(7). The turbulent thermal conductivity is expressed in terms of the eddy 
viscosity by use of the turbulent Prandtl number. Other turbulence models may be 
incorporated into the analysis by suitably modifying or replacing this subroutine. 
If TURB is replaced, the argument list must be identical to the existing list. 

The parameters in the argument list are defined as follov/s. 

KB0DY If KB0DY=1 , forebody/centerbody boundary layer is 

specified. 

If KBj3DY=2, cowl boundary layer is specified. . 

KPLANE If KPLANE=1, initial-value surface is specified. 

If KPLANE=2, solution surface is specified... 

I Circumferential index, of solution point. 

J Radial index of solution point. 

ETX Eddy viscosity for streamwise momentum equation. 

ETZ Eddy viscosity for cross-flow momentum equation, 

ETE Eddy thermal conductivity. 

DEC0DE . This subroutine is used to decode the boundary layer transformed 
variables into physical variables and then to print both the forebody/centerbOdy 
and cowl boundary layer solutions. 

THERMO . This subroutine is used to compute the spatial gradient of the- 
static enthalpy. The assumed thermodynaiTiiG model is that Of a thermally 
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calorically perfect gas. Other thermodynamic models may be incorporated Into 
the analysis by suitably modlfylag or replacing this subroutine. If subroutine 
THERMO Is replaced, the argument list must be Identical to the existing list. 

The parameters In the argument list are defined as follows. 

T temperature 

P pressure 

DTDY first partial spatial derivative of temperature 

DPDY. first partial spatial derivative of pressure 

DHDY first partial spatial derivative of static enthalpy. 

Thermodynamic property Information Is supplied to subroutine THERMO through the 
named common block GASl . 

14. 0VERLAY (4,1) 

LINK41 . This program routine Is the main control routine In OVERLAY (4,1) 
and Is used to control the Integration of the boundary layer differential equa- 
tions.. 

ABLSLN . This subroutine computes the flow properties across the boundary 
layer at a circumferential station that Is on the plane of symmetry. The 

transformed governing differential equations are solved Implicitly using jUie 

matrix solution algorithm In subroutine BTRID. 

AEXT . This subroutine calculates the edge condition parameters for the 
computation of a plane of symmetry boundary layer station. Subroutine AEXT is 
called from subroutine ABLSLN. 

ASPEC . This subroutine determines certain parameters which are used In 
the plane of symmetry boundary layer computation. Subroutine ASPEC Is called 
from subroutine ABLSLN. 

GBLSLN . This subroutine computes the properties across the boundary layer 
at a circumferential station that Is not on a plane of flow symmetry. The trans- 
formed governing equations are solved Implicitly using the matrix solution al- 
gorithm In subroutine BTRID. 

GEXT . This subroutine calculates the edge. condition parameters for the 
computation of a boundary layer station that Is not on the plane of flow sym- 
metry. , Subroutine GEXT Is called from subroutine GBLSLN. 

GSPcC . This subroutine determines certain parameters which are.used In 
the boundary layer computation for stations not oh the plane of flow symmetry. 
Subroutine GSPEC Is called from subroutine GBLSLN. 

BTRID. This subroutine solves a system of simultaneous linear equations 
with a block tridiagonal coefficient matrix. A direct factorization procedure 
employing triangular decomposition Is used to obtain the solution vector. 
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PROD . This subfoutine is used to perform matrix multiplication. Sub- 
routine PRfJD is called from subroutine BTRID. 


15. OVERLAY (4,2) 

LINK42 . This program routine is the main control routine in 0VERLAY (4,2) 
and is used to. determine the boundary layer properties downstream of a shock 

wave-boundary layer interaction region. An integral formulation is used to 

obtain the solution. 

SBLINT . This Subroutine determines the dovmstream boundary layer thick- 
ness and velocity profile exponents in ■ i shock wave-boundary layer interaction 
analysis. . 

INTE6 . This subroutine is used to numerically evaluate certain integral 
parameters involved in the shock v<ave-boundary layer interaction analysis. 


16. 0VERLAY. (4,3) 

LINK43 . This program routine is the main control routine in 0VERLAY (4,3) 
and is used for mesh alteration in the boundary la yer computat ion to account 
for boundary layer growth. 


17. 0VERLAY (4,4) 

LINK44 . This program routine is the main control routine in 0VERLAY (4,4) 
and is used for the internal generation of the cowl boundary layer initial data. 
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SECTION IV 
INPUT PARAMETERS 


Oi^IG?l'-3AL t:;. 
OP POOi< QUAli 1 Y 


1. INTRODUCTION 

The input data required for execution of the computer program are entered 
by both namelist input and formatted read statements. In all cases, the eleven 
namelists LISTl to LISTll are entered. For cases in which the user, selects to 
specify the flow property field on the initial-value plane, that information is 
entered by a formatted read of files ITAPl , ITAP2, and ITAP3. 

In general, only those parameters and data pertinent to the particular prob- 
lem being considered must be entered. Many input parameters have default values 
and do not need.to be specified unless values jather, than the. default values are 
to be entered. 

In this section, each input parameter is defined. Where applicable, both 
the default value and a typical value of the input parameter are given. 


2. TITLE CARD 

The first card of each data deck is a title card on which 72 alphanumeric 
characters (any standard Fortran characters) of identifying information may be 
specified. This card must be the first card of the data deck even if no infor 
mation is listed on it. The format of the card is (12A6). 


3. NAMELIST LISTl 

The parameters entered in namelist LISTl control the overall execution of 
the program. 

KUNIT An integer variable denoting whether English absolute units or SI 

units are to be used in the computation. If KUNIT = 1, English 
absolute units are employed. If KUNIT = 0, SI units are employed. 

A default value of 1 is. specified -for KUNIT. 

KCALL A one-.dimensional integer variable array consisting of three ele- 
ments. Each element of KCALL specifies whether or not a particular 
supersonic flow field integration option is to be performed. If 
KCALL(I) ® 1 (1=1 ,2,3) , then the Corresponding flow field integration 
option is performed. If KCALL(I) = 0, the integration option is not 
performed. The elements of KCALL control the flow field integration 
options and have Specified default values as listed below. 
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XEND 


KTRANS 


KBLAY 


KCALL (I) 

Flow Field Integration Option 

Default Value 

KCALL(l) . 

forebody supersonic flow field . 

1 

KCALL(2) 

incernal supersonic flow field 
with. discrete shock, wave fitting 

1 

KCALL(3) 

internal supersonic flow field 
without discrete shock wave 



fitting 

0 


Specifying KCALL(2) = 1 and KCALL(3) = 1 simultaneously will cause 
an error message to be printed and the program execution to be 
aborted. 


A one-dimensional real variable array consisting of three elements. 
XFND(I) (1=1 ,2,3) denotes the x-position, in either feet. or meters, 
at which the supersonic flow field integration specified by the cor- 
responding element of KCALL is to be terminated. Each element of 
XEND denotes a flow field integration option termination point and 
has a default value as follows. 


XEND(I) 

Termination Point for 

Default Value 

XEND(1 ) 

forebody supersonic flow 
field inte.gration 

2.0 ft 

XEND(2) 

internal supersonic flow 
field integration with 
discrete shock .wave fitting 

3.5 ft 

XEND(3) 

internal supersonic flow 



field integration without 
discrete shock wave fitting 

3.5 ft 


If any element of XEND exceeds the x-position to which the center- 
body geometry is specified, or if XEND(2) or XEND(3) exceeds the 
x-position to which the cowl geometry is specified, appropriate 
error messages are printed and the program execution is aborted. 

Each element of XEND must be positive. It Shoiild be noted. that 
if KCALL(I) = 0, XEND(I) does.not have.to be specified. 

An integer variable denoting whether or not distributed wall bleed 
effects are to be included in the supersonic core flow computation. 
If KTRANS * 1, bleed effects are included in the computation. . If 
KTRANS = 0, bleed effects are not included in the computation. The 
wall bleed distribution is specified by the parameters. entered in. 
namelist LISTS. _A default value of 1 is specified for KTRANS. 

An. integer variable denoting whether dr not the boundary layer com- 
putation is to be performed. If KBLAY = 1, the boundary layer com- 
putation is performed. If KBLAY = 0, the boundary layer computatibn 
is not performed. A default value Of 1 is specified for KBLAY. 
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RCAV6 A positive real variable denoting the estimated average radius, in 

either feet or meters, of the. bow shook wave at XEND(l). The 

specified value of P.CAVG is used in estimating the captured mass 

flow rate at XEND(l). This mass flow rate is used in determining 
if point addition is to be performed in the forebody flow field 
integration. RCAVG must be sped fied Only if the supersonic fore- 
body. flow field integration, option is used [KCALL(l) = 1]. RCAVG 
must be specified even if the forebody integration Option is the 
only integration option used tKCALL(2) = KCALL{3) =0], The cowl 
lip radius may be used for RCAVG. A default value of 0.8 ft is 
specified for RCAVG. 

KVISCY An integer variable denoting whether or not the viscous and thermal 
diffusion terms are to be included in the computation of the fore- 
body supersonic flow field or the internal supersonic flow field. in 
which shock waves are not discretely fitted. If KVISCY * 1, these 
terms are included in the computation. If KVISCY = 0, these terms 
are not included in the computation. A default value of 0 is speci- 
fied for KVISCY. 

KSYM An integer variable denoting the flow symmetry option to be employed 

in the computation. If the boundary layer computation is not employ- 
ed (KBLAY = 0), KSYM can. have the values 0, 1, 2, and 3 corresponding 
to: 

KSYM Flow S.ymmetr.y Option 

0 no planes of symmetry - computed sector is the entire 
solution plane 

1 one plane of symmetry - computed sector is the half- — 

plane bounded by the y-axis and containing the +z- 

axis 

2 two planes of symmetry - computed sector is the 
quadrant bounded by the +y-axis and the +z-axis 

3 axisymmetric flow - computed sector is the single 
circumferential station on the +y-axis 

If the boundary layer computation is employed (KBLAY = 1)» then 
the allowable value of KSYM is 1.. The supersonic flow networks 
for the -flow symmetry, opti ons are illustrated in Figure 3. A 
default value .of 1 is specifi.ed for KSYM. 

KS6L0B An. integer- variable denoting whether or not global correction is to 
be.^jerformed in obtaining the solution for the bow shock wave points. 
If KSGLjDB = 1, global Correction is performed. If KSGLI8B » 0, global 
correction is not performed. Global correction can only be performed 
for the bow shock wave points and not for the internal flow field 
shock wave points. Hence, KSGL0B must be specified only if 
KCALL(I) =1. A default value of 1 is specified for KSGL0B. 
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FIGURE 3. COMPUTATIONAL POINT NETWORKS 
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FIGURE 3. (CONTINUED) 







KPRINT Ar iRteger variable denoting which of tJiree print options is to be 
employed in the program execution for the supersonic flow solution. 
KPRINT can have the following values. 

KPRINT Print Option 

0 print body Solution points. Md shock wave solution 
points only 

1 print all solution points „ 

2 print all solution points and shock wave parameters 

(incident normal Mach number and shock wave surface 
normal unit vector components) 

A default value of 1 is specified for KPRINT. 

IPRSTP A positive Integer variable denoting the plane number at which the 
program execution is to be terminated. This input parameter is 
typically employed when the execution is intended to be restarted 
at plane number (IPRSTP + 1). Specifying IPRSTP < 0 has no effect 
on the execution, of the program. A default value of 0 is specified 
for IPRSTP. 

KSTART An integer variable used in controlling the file operations involved 
in the restarting of the program (if desired). KSTART controls the 
storage and retrieval of the restart file as follows. 

KSTART Control of Restart File 

0 no file operations 

1 write restart file on TAPE 9 

2 read information for program restart from TAPE 9, 

and write ensuing solution planes, on this tape 

A default value of 0 is specified for KSTART. TAPE 9 is linked to 
the dummy file RESTRT in the PROGRAM card. 


4. NAMELIST LIST2 

The parameters entered in namelist LIST2 specify the free-stream condi- 
tions, the inlet orientation, and the parameters which control the internal 
generatioti af the supersonic flow property field on the initial-value plane. 

A positive real variable denoting the free-stream Mach number. The 
specified value of MFS must be greater than 1.0. A default value of 
3.0 is specified for MFS. 

A positive real. variable denoting the free-stream absolute pressure. 


MFS 

PFS 
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R(9FS 


PITCH 


YAW 


XI, 


KIVS 


KC0N 


In 6lth6r (1bf/ft^) or (N/m^)# A default valud of 242,2 (Ibf/ft*^) 
is specified for PFS (this value is the pressure of the standard 
atmosphere at an altitude of 50,000 ft). 

A positive real variable denoting the free-stream density, in either 
(slug/ft3) or (kg/m3). A default value of 0.0003622 (slug/ft3) is 
specified for R0FS (this value is the density of the standard atmos- 
phere at an altitude of 50,000 ft). 

A real variable denoting the angle, in degrees, subtended by the free- 
stream velocity vector and the projection of the free-stream velocity 
vector on the (x,z)-plane, as illustrated in Figure 4. A default 
value of 1.0 degree is specified for PITCH. If KBLAY-1 is specified 
in namelist LISTl, then PITCHS^O.O. 


A real variable denoting the angle, in degrees, subtended by the x- 
axis and the projection of the free-stream velocity vector on the 
(x,z)-plane, as illustrated in Figure 4. A default value of 0.0 
degrees is specified for YAW.. 

A positive real variable denoting the axial, (x) position, in either 
feet or meters, of the supersonic flow initial-value plane. If the 
forebody flow field integration option is specified [KCALL(l) =1J, 
xr must be specified at the beginning of the forebody flow field 
computational regime (see Figure 1). If only the internal flow 
field integration option is specified [KCALL(l) = 0, KCALL(2) = 1 
or KCALL(3) = 1], XI must be specified at the axial station of the 
cowl lip (see Figure 1). XI must not fall outside of the range of 
axial stations .for which the centerbody geometry is specified. Also, 
XI must not be greater than the axial station up to which the cowl 
geometry is specified. A default vaule of 1.0 ft is specified for 
XI. 


An integer variable denoting whether the supersonic flow field on 
the initial-value plane is to be generated internally or read in. If 
KIVS = 1, the initial-value plane data are computed internally. If 
KIVS = 0, the initial-value plane data must be supplied by the user 
through a formatted read of file ITAPl . The formatted read input is 
described at the end of this section. The internally generated initial 
value plane option is applicable only to cases where the forebody is 
conical up to the axial station where the initial-value plane Js lo- 
cated. Specifying KIVS =1 requires that YAW = 0.0. A default value 
of 1 is specified for KIVS. 


An integer variable denoting whether or not the bow shock wave is 
conical, KC0N must be Specified only if the forebody flow field in- 
tegration option is employed [KCALL(l) =1]. If KC0N = 1 , the bow 
shock wave is conical. In this case, the angle at each initial-value 
plane bow shock wave point that. is Subtended by the shock wave and 
the x-axis in the meridional plane defined by the shock wave point ^ 
is computed internally. If KG0N = 0, the bow shock wave is not coni- 
cal, and the Shock wave angles must be supplied by the user through a 
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formatted read of file ITAPl. The formatted read Input Is described 
at the end of this section. If the initial-value plane is generated 
Internally (KIVS * 1),, a conical bow. shock wave Is assumed (KC0N = 1). 

A default value of 1 Is specified for KC0N. 

KSUPER An Integer variable denoting the method to be used In Internally gener- 
ating the supersonic flow Initial data. KSUPER must be specified only 
If the Initial data are to be generated Internally (KIVS =1). If 
KSUPEPv = 1, the supersonic flow Initial data are computed using the. 
approxlmate Taylor-Maccoll algorithm described in Section II. If 
KSUPER = 2, the supersonic flow initial data are computed using the 
Jones algorithm described In Section II. If KSUPER = 1, KSYM Can 
have the values 1, 2, or 3, where if KSYM = 3 the angle of attack 
must be 0.0 (PITCH = 0.0). If KSUPER = 2, KSYM must be equal to 1 . . 
and PITCH f 0.0. A default value of 2 Is specified for KSUPER. 

ITAPl A positive Integer variable denoting the tape number from which the 

user supplied supersonic flow Initial data are to be. entered by a 
formatted read. ITAPl must be specified only If the Initial-value 
plane flow property field Is to be supplied by the user (KIVS * 0). 

The default value assigned to ITAPl is 5 (the Input file). The 
user may specify ITAPE = 10, In which case the supersonic Initial 
data are read from TAPE 10. TAPE 10 Is linked to the dummy file 
IVSl in the PR0GRAM card. 


5. NAflELIST LISTS 

The parameters entered in namelist LISTS specify the number of circumfer- 
ential and radial stations used in the supersonic flow computational point net- 
work. 

IST0P A positive integer variable denoting the number of circumferential 
stations used in the supersonic solution computed flow field sector. 

The value specified for IST0P must correspond to the flow symmetry 

option specified by KSYM (see namelist LISTl) as follows (see Figure 3). 

KSYM Allowable Value(s) of IST0P 

0 5 < IST0P < 30 

1 4 < IST0P < 16 

2 3 < IST0P < 8 

3 IST0P = 1 

A default value of 15 is. specified for IST0P (this recommended value 
corresponds to KSYM = 1). If the supersonic initial data are to be 
generated internally using the Jones algorithm (KIVS = 1, KSUPER =2), 
then IST0P must be greater than 1. 
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JMAXI 


JINLET 


JLIMIT 


Integer variable denoting, the number of radial stations 
on the supersonic flow Initial-value plane. Note that the Inltlal- 

stlrtPri format Is the same whether the computation Is being 

ni. axIaVstatlon upstream of the forebody computational 

flow regime or at the axial station of the cowl lip. The specified 
value for JMAXI must be at least 3 and nO greater than 15. The de- 
fault and recommended value for JMAXI is 11. 

JMAXI, the outermost radial station corresponds 
lldut shock wave points, and the remaining (JMAXI - 1) 

raaial stations correspond to the streamline points. 

denoting the number of radial stations on 
each solution plane in the supersonic internal flow field Inte- 

specified only If the internal flow field 
Is to be computed. The specified value of JINLET must be at least 

o/jINLET^^fir value and the default value 

r^r inteqratlon option In which shock waves 

are not discretely fitted [KCALL(2) =0, KCALL (3) = 1] JINLET 

station*^ For streamline points at each circumferential 

Station. For the internal flow field integration option in which 

“5“'^ ““®! /te discretely fitted [KCAU(2) = 1, KCALL(3) =0], the 

?o f ?* f'-tuinferentlal station Is equal 

to (JINLET - 2). The remianing two storage locations are assigned 
to the upstream and downstream shock wave ..solution points. ^ 

n variable array consisting of two elements. 

JLIMIT must be specified only if the forebody flow field integration 

[KCALL(l) = 1]. The elements of JLIMIT are used 
in controlling the number of Interior field points which are added 

hmi nr?S??nn 2^ The fl^st eUmenl of 

J IMn [JLIMIT(l)] denotes the allowable maximum number of radial 

® solution plane in the forebody flow field computation 
when the mass flow rate across that plane is less than a specified 
fraction [denoted by CRIT(7) In namelist LIST10]of the estimated ° 
flo^fip?d^?Jf«^^ XEND(1), the axial location at which the forebody 

terminated. The second element of 
'^1:^5.. LJLpIT(2)] denotes the allowable maximum number, of radial 
stations when the mass flow rate exceeds the specified fraction of 
the estimated mass flow rate at XEND(l). Each element of JLIMIT 

^ '10 greater than 15. 

JLIMIT(l) should be less than or equal to JLIMIT(2). A default 
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6. NAMELIST LIST4 


The parameters entered in namelist LIST4 specify the thermodynamic model 
and the molecular transport, properties. 

R A positive real variable denoting the gas constant, in either 

(ft-lbf)/(slug-R) or (J/k.g-K). A default value of 1716Jill6 
{ft-lbf)/(slug-R) is Specified for R. 

GAMMA A positive real variable denoting the specific heat ratio. A 

default value of 1.4 is specified for GAMMA. 

If the boundary layer computation is not employed (KBLAY = 0) and if the 
influence of molecular transport, is not to be included in the supersonic flow 
computation (KVISCY * 0), then no other parameters need be specified in namelist 
LIST4. If these options are employed in the computation, then the parameters 
presented in the following discussion must be specified. 

The molecular dynamic viscosity is represented in the computer program by 
the Sutherland formula (6). 



In equation (1), y is the dynamic viscosity at the absolute temperature T, and 
B is a constant. For air, B has the value of 198.6 R (110 K). The parameter 
Uq is the. viscosity at the reference temperature Tq. The constants Po» Tq* 

B must be specified in the program input by entering the following three para- 
meters. 

VIS0 A positive real variable denoting the reference viscosity uq ii^ 

equation (1). The units of VIS0 are either (lbf-sec/ft2) or 
(N-s/m2). A default value of 3.5 x 10*7 (lbf-sec/ft2) is specified 
for VIS0 (this value is the dynamic viscosity of air at 492.0 R). 

T0 A positive real variable denoting the reference absolute. temperature 

Tq in equation (1.). The units of T0 are either R or K. The specified . 
value, of T0. must correspond to the specified value of VIS0. A de- 
fault value of 492.0 R is specified for T0. 

6 A positive real variable denoting the constant B in equation (1). 

The units for B.are either R or K. A default value of 198.6 R is 
specified for B. 
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The molecular thermal conductivity is represented in the computer program 
by 

K = CpU/Pr 

where k denotes the thermal conductivity, Cp is the constant pressure specific 
heat, and Pr is the laminar Prandtl number which is assumed constant in the 
analysis. The specific heat Cp is calculated internally in the program, whereas 
the Prandtl number Pr is user s|iecified by entering, the following parameter. 

PR A positive real variable denoting the laminar Prandtl number. A 

default value of 0.71 is specified for PR. 


7. NAMELIST LISTS 

The input parameters entered in namelist LISTS specify the contours of the 
centerbody and the cowl. It is assumed that both the centerbody and the cowl 
are axisymmetric. The x-coordinate axis is the longitudinal axis of both the 
centerbody and the cowl (see Figure 1). The forebody tip must be located at 
X = 0.0. The axta.1 station of the cowl lip must be specified at x > 0.0. 


For the purpose of geometry description, the axial (x) domain is divided 
Into a number of intervals, as illustrated in Figure S. The '^^’’^ber of axial 
stations at which the centerbody geometry is specified is denoted by NCENT. The 
number of intervals on the centerbody is equal to (NCENT - 1). The number 
axial stations at which the cowl geometry is specified is denoted by NC0WL. The 
number of intervals on the cowl is equal to (NC|3WL - 1). 


In any interval, the centerbody or cowl radius may be specified either 
by tabular input, or by supplying the coefficeints in a cubic polynomial written 
as a function of x. For the tabular input option, the body radius rU) at axial 
position X in the i th interval (xi^^ x < xt+i) is found by linear interpolation 
between point (xi,riT and point (xi+l,ri+i). The local slope ^°r. 

this interval in a given meridional plane is then given by the slope of the line 
segment joining these two points. Alternatively, employing the cubic polynomial 


r(x) = a^ + b^(x - x^) + c^(x - x^)^ + d^(x - x^)^ (x^ ^ x < (3) 

requires that the curve fit coefficients ai,. bi, Ci , and di be supplied by the 
user. Since equation (3) is a cubic, slope and curvature can be matcned at the 
junction point between two adjacent intervals employing this formulation. An 
option exists to employ the following cubic polynomial instead of equation 

r(x) s + b^x + c^x^ + d^x^ (x^ < x < x^^^ ) (4) 
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1st INTERVAL End INTERVAL (NC0WL-l)th INTERVAL 



FIGURE 5. BOUNDARY CONTOUR DESCRIPTION 


KBASE 


f or.equatlon (4), the coefflefents s,. b,, c. and 

Tile a^a? ' ' ) anO/oT (NC0wi - 1) tJtehvaU. 

always be specified for NCENT and (if the internal 
U°used^ ^ computed) NC0WL axial stations, no mattec_wlilch_-formulation 

rAnfoiLn^^ forebody is conical ahead of a certain axial station, the forebody/ 

ind ?^A ?LaT?^^ interval) may be specified by enter- 

ing the cone half-angle directly rather than by supplying the curve fit coeffi- 
cients or entering the body radius by tabular input. 

user Il!d HopfJl? description option for a given interval is. specified by the 
fo^hnrfw^AnL 1 ° intervals.. For instance, the 

IaHa?inn ® conical tip. then a quadratic or cubic 

cSbiflaHsHl ® linear variation with x, then again a quadratic or 

describe this contour would be to input the 
cone half-angle for the first interval, the cubic curve fit coefficients for 
the second interval, the body radii at the ends of the third interval, and the 
coefficients for the fourth interval. Alternatively, the appro- 
coefficients could be supplied for each of the four in- 
tervals. It should be noted that in the input of the geometry data, radius 
slope, and curvature should be made compatible between adjacent intervals The 
boundary contours are specified by entering the following pa Seters 

An integer-variable denoting whether equation (3) or equation (4) is 
0 e einployBd when at least one interval of the centerbody or cowl 
geometry is specified by either of these two cubic equations. If 
neither equation (3) nor equation (4) is used in the geometry de- 
scription, then KBASE does not have to be specified. If a cubic 
equation is to be used, then entering KBASE = 0 specifies that 

employed. If KBASE = 1, then equation (4) will 
° specified for KBASE. The specified 
value of KBASE applies to both the forebody/centerbody contour and 
the cowl contour. 

A positive integer variable denoting the number of axial stations 
used in specifying the forebody/centerbody geometry. The number 
of intervals for the forebody/centerbody is equal to (NCENT - 1). 

The specified value for NCENT must be at least 2 and no greater 
than 25. A default value of 2 is specified for NCENT. 

A one-dimensional -integer variable array (dimensioned at 25) consist- 
ing of. (NCENT -1) elements. Each element of KDCENT specifies the 
forebody/centerbody geometry description option to be .used in the 
corresponding^interval. Specifying KDCENT( I) = 1 for the i th inter- 
val selects the option in which the forebody/centerbody radius is 

or equation (4) (depending on the val^e of 
KBAp), Specifying KDCENT(l) » 2 for the i ^h interval Selects the 
option in which the forebody/centerbody radius is specified by tabu- 
lar input. If the forebody tip is conical, the geometry in the first 
forebody/centerbody interval may be described by entering the cone 


NCENT 


KDCENT 
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half-angle directly and specifying KDCENT(l) “ 3. The specified 
value for KDCENT(I) must be 1 or 2, except for the first interval 
(I » 1) in which Case values of 1, 2, or 3 may be specified. Speci- 
fying other A^alues than those allowed causes the program execution 
to be aborted. The default value for KDCENT(l).is 3, while all _ 
other elements of KDCENT have.no default values specified. 

XCENT A one-dimensional real variable array (dimensioned at 25) consisting 

of NCENT elements. Each element of XCENT denotes the axial (x) posi- 
tion, in either feet or meters, of the beginning of a forebody/center- 
body interval [XCENT(NCENT) denotes, the axial position of the end of 
the last interval]. The elements of XCENT must be nonnegative and 
monotonically incresing. The default values for XCENT(l) and XCENT(2) 
are 1.0 ft and 3.5 ft, respectively. The remaining elements of XCENT 
do not have default values specified. 

RCENT A one-dimensional real, variable array (dimensioned at 25) consisting 
of up to NCENT elements. Each element of RCENT specifies the fore- 
body /centerbody radius, in either feet or meters, at the axial loca^ 
tion specified by the corresponding element of XCENT. If KDCENT(I) = 2, 
then RCENT(I) and RCENT(I + 1) must. be specified. If KDCENT(I) = 1 or 
3, then RCENT(I) and RCENT(I +1) do not have to be specified. Each 
element of RCENT must be nonnegative. No default values are specified 
for the. elements of RCENT. 

AGENT One-dimensional real variable arrays (each array is. dimensioned at 

BCENT 25), where each array consists_of up to (NCENT - 1) elements. These 

CCENT arrays are used in conjunction with equation (3) or equation (4) for 

DCENT specification of the forebody/centerbody geometry. The elements of 

AGENT,, BCENT, CCENT, and. DCENT specify the coefficients ai , bi, Ci , 
and d|, respectively, in equation (3) or equation (4). If KDCENT(I) = 

1, then ACENT(I), BCENT(I), CCENT(I), and DCENT(I) must be specified 
for the ith interval. If KDCENT(I) = 2 or 3, then ACENT(I), BCENT(I), 
CCENT(I), and DCENT(I) do not have to be specified for that Interval. 

The units for the elements of AGENT are either feet or. meters. The. 
elements of BCENT are dimensionless. The units for the. elements of 
CCENT are either (feet)“l- or (meters)“i. The units for the elements 
of DCENT are (feet)“2 or (meters)“2. No .default values are specified 
for the elements of AGENT, BCENT, CCENT, and DCENT. 

C0NE A real variable denoting the cone half-angle, in degrees, of the 

forebody tip if it is conical. If KDCENT(l) * 3, then C0NE must.be 
specified. If KDCENT.(1) = 1 or 2, then C0NE does not have to be 
specified. A default value of 10.0 degrees is specified for C0NE. 

If only the externdl.-flow field about the forebody is. to be Computed 
[KCALL(l) = 1, KCALL(2) » 0,.and KCALL(3) * 0], no further parameters must be 
specified in namelist LISTS. If the flow field in the annulus is to be deter- 
mined, the_follOwing parameters must be entered. 


NCfSWL 


KDC0WL 


XC0WL 


RC0WL 


AC0UL 

BC0WL 

CC0WL 

OC0WL 


DXTRAN 


A positive Integer' variable denoting the number of axial stations 
used In specifying the cowl geometry. The number of Intervals foi* 
the cowl Is equal to. (NC0W1L - 1). The specified value for NC0WL 
must be at least 2 and no greater than 25. No default value Is 
specified for NC0WL. 

A one-dimensional integer variable array (dimensioned at 25) consist- 
ing of (NC0WL - 1) elements. Each element of KDC0WL specifies the 
cowl geometry description option to be used for the corresponding 
Interval. Specifying KDC0WL(I) = 1 for the 1 interval selects 
the option in which the cowl radius is described by equation (3) or 
equation (4) (depending on the value of XBASE). Specifying KDC0WL(I) = 
2 for the i ^ Interval selects the option in which the cowl radius, is 
specified by tabular input. The specified value for KDC0WL(I) must be 
either 1 or 2. Specifying other values than those allowed causes the 
program execution to be aborted. The elements of KDC0WL do not have 
default values specified. 

A one-dimensional real variable array (dimensioned at 25) consisting 
of NC0WL elements. Each element of XC0WL specifies the axial (x) 
position, In either feet or meters, of the beginning of a cowl inter- 
val [XC0WL(NC0WL) denotes the axial position of the end of the last 
interval]. Each element of XC0WL must be nonnegative and monotonically 
increasing.. The elements of XC0WL do not have default values specified. 

A one-dimensional real variable array (dimensioned at 25) consisting 
of up to_NC0WL elements. Each element of RC0WL specifies the cowl 
radius. In either- feet or -meters, at the axial location specified by 
the. corresponding element of XC0WL. If KDCOWL(I) = 2, then. RC0WL(I) 
and RC0UL(I + 1 ) must be specified. If KDC0WL(I) = 1, then RC0WL(I) 
and RC0WL(1 + 1 ) do not have to be specified. Each element of RC0WL 
must be nonnegative. The elements of RC0WL do not have default values 
speci fied. 

One-dimensional real variable arrays (each array is dimensioned at 
25), where each array consists of up to (NC0WL - 1) elements. These 
arrays are, used in conjunction with equation (3) or equation (4) for 
specification of the cowl geometry. The elements of AC0WL, BC0WL, 

CC0WL, and DC0WL specify the coefficients a^ , b<, cs, and di, respec- 
tively, In equation (3) or equation (4). If KD60WL(I) = 1, then 
AC0WL(I), BC0WL(I), CC0WL(I) , and. DC0WL(I) must.be specified for the 
Itb Interval. If KDC0WL(I) = 2, then AC0WL(I), BC0WL(I) ,. CC0WL(I) , 
and DC0WL(I) do not have to be specified for that Interval. The 
units of the elements of AC0WL are either feet Or meters. The ele- 
ments of BC0WL are dimensionless. The units for the elements of 
CC0WL are either (feet)-l or (meters)-i. The units for the elements 
Of DC0WL are either (feet)~2 or (meters)~2. No default values are 
sped fled, for the elements of AC0WL, BC0WI., CC0WL, and DC0WL. 

A real variable denoting the cehterbody translation from the design 
point position, or, equivalently, the amount the cOwl has been trans- 
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lated with respect to the centerbody. The units of DXTRAN are- 
elther feet or meters. Translation occurs solely in the x-direction. 
Moreover, the origin the -coordinate system is maintained at the 
forebody tip. when translation Occurs. A positive value for DXTRAN 
corresponds to a forward centerbody translation or a rearward cowl 
translation. A default value of 0.0 is specified for DXTRAN. 


8. NAMELIST LIST6 

The parameters entered in namelist LIST6 control the internal generation of 
the boundary layer initial data and specify the boundary layer computational mesh. 
The parameters entered in this namelist need only be specified if the boundary 
layer computation is to be performed [KBLAY » 1 specified in namelist LISTl]. 

The following six parameters are used for specification of the forebody/ 
centerbody boundary layer initial data. 

KBLIDA An integer variable denoting whether or not the forebody/centerbody 
boundary layer initial data are to be internally generated by the 
Adams, finite .difference algorithm described in Section II. The Adams 
algorithm is applicable only to cases in which the forebody is conical 
ahead of the station where the boundary layer data are to be generated. 
If KBLIDA = 1, the initial data are internally generated. If KBLIDA = 

0, the initial data are entered by a formatted read of file ITAP2. The 
formatted read is described at the end of this section. A default value 
of 1 is specified for KBLIDA. 


If KBLIDA = 1, the following four parameters must be entered. . 

KTURB An integer variable denoting if the Adams algorithm is to generate 

laminar flow or turbulent flow initial data for the forebody/center- 
body boundary layer. If KTURB = 0, laminar flow initial data are 
generated. If KTuRB - 1, turbulent flow initial data are generated. 

A default value of 0 is specified fOr KTURB. 

CLEN6H A positive real-variable denoting the length of the conical section. 

of the forebody in either feet or meters. A default value of 4.0 
ft is specified. for.CLENGH. 

The following two parameters define the Computational mesh used in the Adams 
algorithm solution [see Reference (5) for further discussion.]. In general, the 
pCOgram is executed by retaining these input parameters at their default values. 

ADY A positive real variable denoting the first transformed mesh length 

in the Surface normal direction that is used in the Adams algorithm 
Solution — The default and. recommended value of AD.Y is 0.010. 

ARATI0 A positive real variable denoting the ratio of successive normal 
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mesh steps used in the Adams algorithm solution. The default and 
recommended value of ARATI0 is 1.0630. 


If KBLIDA = 0 is speci fied, in this namelist, then the following parameter 
must.be entered. 

ITAP2 A positive integer variable denoting the tape i.'imber from which the 

forebody/centerbody boundary layer initial data are to be entered 
by a formatted read. The default value assigned to ITAP2 is 5 (the 
input file). The user may specify ITAP2 = 11, in which case the 
forebody/centerbody boundary layer initial data are read from TAPEll. 
TAPEll is linked to the dummy file IVS2 in the PROGRAM card. 

The following five parameters are used for specification of the cowl boun- 
dary layer initial data and must be entered only if the cowl boundary layer ir 
to be calculated [KCALL(2) = 1 or KCALL{3) = 1 and KBLAY = 1 specified in nar.e- 
list LISTl]. 

KBLID3 An integer variable denoting whether or not the cowl boundary lay'ir 
initial data are to be internally generated by an approximate tech- 
nique described in this section. If KBLIDB = 1, the initial data 
are internally generated. If KBLIDB = 0, the initial data are 
entered by a formatted read of file ITAP3. The formatted read is 
described at the end of this section. A default value of 1 is 
specified for KBLIDB. 

The internally generated cowl boundary layer initial data are obtained using 
an approximate analysis briefly described below and presented in g>"eater detail 
in Reference (1). The initial step in. the analysis is to define the boundary 
layer thickness at each circumferential station in the computed sector. For 
l.aminar flow, the lacal boundary layer thickness 6 is approximated by 


6 = 5.0 



0 u 
*^e e 


1/2 


(5) 


whereas for turbuJent flow, 6 is approximated by 


6 ^ 0.37X 


PgUeX 


-1/5 


( 6 ) 


where pe* Ue* *^e denote the local boundary layer edge density, streamwise 
velocity, and laminar dynamic viscosity, respectively, and”? is the streamwise 
eurvi* inear coordinate on the cowl (see Figure 5). With the boundary layer 
thickness determined, the streamwise and croSS-flow velocity profiles for laminar 
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ORIGINAL RAGK U-S 
OF POOR QUALITY 


flow are defined by 
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for turbulent flow, the following power laws are employed 
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In the ebo.e equetlons. u end i denote 

la the Tall, a iJd Tare constants^and the subscript e denotes boundary 

Uver edge conditions. The total enthalpy H distribution is given by assuming 
a quadratic variation with distance across the boundary layer and then imposing 
the appropriate boundary conditions. This yields 


H - + 2(Hg - 




+ (H„ - H^) 




(11) 


for cases in which the wall temperature .is specified, and 


H • (Hj - 
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for «,ds in Which ‘!|| ?oUl!d“?h’alM ird'»ta?’en{"h.lpy 
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profiles allow the density distribution to be determined if the pressure is 
assumed to be constant along a.^iven boundary layer normal. 

The following three parameters must be entered if KBLIDB =5 1 and are used 
for specification of the cowl boundary .layer initial data. 

KTURBB An integer variable denoting whether laminar flow or turbulent flow 
cowl boundary layer initial data are to be generated internally by 
the approximate analysis presented above. If KTURBB * 0, laminar, 
flow initial data are generated. If KTURBB = 1., turbulent flow 
initial data are generated. A default value of 0 is specified for 
KTURBB . 

ALPCWL A positive real variable denoting the exponent a in the turbulent 

power law streamwise velocity profile given by equation (9). ALPCWL 
must be entered only if KTURBB - 1 is specified in this namelist. A 
default value of 1. 0/7.0 = 0.142857 is specified for ALPCWL. This 
value of ALPCWL ocrresponds to the assumptions made in deriving 
equation (6). 

BETCWL A positive real variable denoting the exponent 3 in the turbulent 
power law cross-flow velocity profile given by equation (10). 

BETCWL must be entered only if KTURBB = 1 is specified in this 
namelist. A default value of 1. 0/7.0 = 0.142857 is specified for 
BETCWL. This value of BETCWL ocrresponds to the assumptions made 
in deriving equation (6) . . 

If KBLIDB = 0 is specified in this namelist, then the following parameter 
must be entered. ^ 

ITAP3 A positive integer variable denoting the tape number from which the 

cowl boundary layer initial data are to be entered by a formatted 
read. The default value assigned to ITAP3 is 5 (the input file). 

The- user may specify ITAP3 = 12, in which case the cowl boundary 
layer initial data are read from TAPE12. TAPE12 is linked to the 
dummy file I VS 3 in the PR0GRAM card. 


The following^two parameters designate the values of the streamwise curvi- 
linear coordinate x at. which the boundary layer initial data are specified. 

See Figure 5 for definition of the x coordinates for both the forebody/center- 
body and cowl. Note that for the forebody/centerbody, x = 0.0 is located at 
the forebody tip, while for the cowl, x » 0.0, is located at the cowl lip. 

XASTRT A positive real variable denoting the value of the streamwise curvi- 
linear coordinate X (see Figure 5) at which the forebody/centerbody 
boundary layer initial data are specified. XASTRT must be entered 
only if KBLIDA = 0 is Specified in this namelist. XASTRT must cor- 
respond to the value of XI specified in namelist LIST2. No default 
value is specified for XASTRT. 


53 


XBSTRT 


°,l ‘'’^ curvl- 

initial data arp specified ^ the cowl boundary layer 

Ms sp«if1ar?n thf! *f ^BLIM . 

location of the first solut1on*olan*f IJIidf 

value IS specified for XBSTRT ^ s de the_annulus. No default 

ti onarLIh! '°The\M„"da?7"a“e?%ompot1?l‘'on‘aI 

Figure 6. ^ ^ computational point networks are illustrated in 

stationrL^brused\^rthe\omputa^ circumferential 

body and cowl boundary layers forebody/center- 

least 5. but no greateT tha^ 6^ at 

value of MBLAY for KSYM = 1 is * 15 ^^^ default value and recommended 

lowing two parameters arrentered^o^deslgnat^th^^ 

the supplied data. oesignate the number of radial points in 

the user-Lppliirfo?eb?dy/ce^^ in 

NINA must be entered only if KBLn?^- boundary layer initial data. 

^ru'L'r-s^ppl?Id'cow[^^^^^^ points in . 

but no greater than 20 A defau?t"^of^?n‘ • ^ 

uiari CU. A default of 20 IS specified for NINB. 

the boundary^ lay3r computSo^^^ Points used in 

KINA and NINB since the oser-sipplied dM^Jf/SeTte^poUted? f™"' 

lploye”fn’?he’f^re'bSdy/«„tl?SJ^^ stations 

ulpeli^JeJ-ffr ^a"‘ 


NB 


enpIoyed*in'the^cowl*bouhdarv* 1 aeI''^ number of radial stations 
fied ;„Iy if the cowl boJnda^v ?aeer‘?^ ““ I”® *P«'- 

at least 5 but no greater than 2(/ A 

fied for NB. y^«<*^or tnan 20 . A default value of 20 is.speci- 


54 


OF pQOn i 


BOUNDARY 

LAYER 

STATIONS^y 


<f>~w 


IB=MBLAY 


COWL- 


CENTERBODY 


lB=i. 




z,vy 


(a) CIRCUMFERENTIAL NETWORK 


JB=NA. NB, 
NINA, NINB 



(b) RADIAL NETWORK 

FIGURE 6. BOUNDARY LAYER COMPUTATIONAL 
POINT NETWORKS 



The following parameter specifies the mesh distortion in the boundary 
layer normal direction. 

FACTOR A one-dimensional real variable array consisting of two elements.- 

The elements Of FACT0R denote the ratio -of two successive step sizes 
in the boundary layer tiormal direction, the normal step size being 
expressed ia terms of transformed variables [see Reference (1)]. 

FACTOR (1) and FACTOR (2) denote the successive normal step size 
ratios for the forebody/centerbody and cowl boundary layers, respec- 
tively. FACTOR (2) must be entered only if the cowl- boundary layer 
is to be calculated. Specifying FACTOR(I) - 1.0 produces a mesh 
that is uniform in the boundary layer normal direction, whereas speci- 
fying FACTOR(I) > 1.0 produces a mesh which is more closely spaced 
near the wall. Default values of 1.0 are specified for both elements 
of FACTOR. 


9. NAMELIST LIST? 

The parameters entered in namelist LIST? specify the wall temperature boun- 
dary condition to be used in the boundary layer computation. The parameters 
entered in this namelist need only be specified if the boundary layer computation 
is to be performed (KBLAY = 1 specified in namelist LlSTl). 

KTYPE An integer parameter denoting whether the wall temperature or wall 

temperature normal derivative is to be used as a boundary condition 
in both the forebody/centerbody and cowl boundary layer solutions. 

If KTYPE * 1, the wall temperature is used as a boundary condition. 

If KTYPE = 2, the normal derivative of the temperature at the wall 
is used as a boundary condition. If the forebody/centerbody boun- 
dary layer initial data are generated internally (KBLIDA = 1 speci- 
fied in namelist LISTS), then the wall temperature boundary condi- 
tion option must be specified (KTYPE =1). A default value of 1 
is specified for KTYPE. 

The wall temperature or its normal derivative may be specified as being 
either a constant or as having functional dependence upon axial position. For 
the case of axial position dependency, the wall temperature or its derivative 
is specified by tabular input. The temperature boundary condition .on the fore- 
body/centerbody is specified by entering the following parameters. 

KWLTA An integer variable denoting whether a constant or variable tempera- 

ture or temperature derivative boundary condition is to be used .for. 
the forebody/centerbody boundary layer computation. If KWLTA * 1, 
a constant boundary condition is specified. If KWLTA =2, a variable 
boundary condition is specified. A default value of 1 is specified 
for KWLTA. 

If a constant temperature or temperature derivative boundary condition is __ 
to be specified, then one of th e two following parameters is entered. 
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TC0NSTA A positive real variable denotiog the constant wall temperature, 
in either. R or K, for-the forebody/centerbody. TC0NSTA must be 
entered only if KTYPE = 1 and KWLTA « 1 are Specified in this name- 
list. The free-streiim stagnation temperature is recommended for 
TC0NSTA. A default value of 500.0 R i^ specified for TC0NSTA. 

DTDYCA A real variable denoting the. constant wall temperature normal 
derivative, in either (R/ft) or (K/ni), for the forebody/center- 
body. DTDYCA must be entered only if KTYPE = 2 and KWLTA = 1 
are specified in this namelist. A default value of 0.0 is speci- 
fied for DTDYCA. 

If a varibale temeprature or temperature derivative boundary condition is 

to be specified, then the following parameters are entered. 

NTABA A positive integer .variable denoting the nun4>er of tabular data 

points used in the temperature boundary condition specification on 
the forebody/centerbody. NTABA must be entered only if KWLTA = 2 
is specified in this namelist. NTABA must be at least 3 but no 
greater than 25. No default value is specified far NTABA. 

XWLA A one-dimensional real variable array (dimensioned at 25) consisting 

of NTABA elements. The XWLA array must be entered only if KWLTA = 2 
is specified in this namelist. Each element of XWLA specifies the . 
axial (x) position, in either feet or meters,, at which the wall 
temperature or its normal derivative is to be specified for the 
forebody/centerbody. Each element of XWLA must be nonnegative, 
monotonically increasing, and within bounds of the specified geom- 
etry. No default values are specified for the elements of XWLA. 

TWLA A one-dimensional real variable array (dimensioned at 25) consisting 

of NTABA elements. The TWLA array must be entered only if KTYPE = 1 
and KWLTA = 2 are Specified in this namelist. Each element of TWLA 
specifies the forebody/centerbody wall temperature, in either R or K, 
at the axial position Specified by the corresponding element XWLA. 
Each element of TWLA must be positive. No default values are speci- 
fied, for the elements of TWLA. 

DTDYWA A one-dimensional real variable array (dimensioned at 25) consisting 
of NTABA elements. The DTDYWA array must be entered only if KTYPE - 
2 and KWLTA = 2 are Specified in this namelist. Each element of- 
DTDYWA specifies the forebody/centerbody wall temperature normal 
derivative,. in either (R/ft) or (K/m)., at the axial position speci-. 
fied by the corresponding element of XWLA. No default values are 
specified for the elements of. DTDYWA. 

If only the forebody flow field and boundary layer are to be calculated 

[KCALL(l) 1 , kCALL( 2) = KCALL(3) = 0^ and KBLAY = 1 Specified in namelist 

LiSTl],, then no other input parameters have to be specified in this namelist. 

If the internal, flow field and boundary layers are to be computed CkCALL(2) = 

1 or KCALL(3) = 1, and KBLAY * 1 specified in namelist LISTIJ, then the follow- 
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Ing parameters must be entered to specify the wall temperature boundary condition 

for the cowl boundary layer_solution. - 

KWLTB An integer variable denoting whether a constant or variable tempera- 

ture or temperature derivative boundary condition is to be used for. 
the cowl boundary layer computation. If KWLTB - 1, a constant boun-. 
dary condition is specified. If KWLTB = 2, a variable boundary con- 
dition is specified. A default value of l .is specified for KWLTB. 

If a constant temperature or temperature derivative boundary condition is 

to be specified, then one of the two following parameters is entered. 

TCjSNSTB A positive real variable denoting the constant wall temperature, in 
either R or K, for the cowl. TC0NSTB must be entered only if 
KTYPE = 1 and KWLTB * 1 are specified in this, namelist. The free- 
stream stagnation temperature is recommended for TC0NSTB. A de- 
fault value of 500. OR is specified for TC0NSTB. 

DTDYCB A real variable denoting the constant wall temperature nor- 
mal derivative, in either (R/ft) or (K/m), for the cowl. DTDYCB 
must be entered only if KTYPE =• 2 and KWLTB = 1 are specified in 
this namelist. A default value or 0.0 is specified for DTDYCB. 

If. a variable temperature or temperature derivative boundary condition is 

specified, then the following parameters are entered. 

NTABB A positive integer variable denoting the number of tabular data 

points used in the temperature boundary condition specification .on .. 
the cowl. NTABB must be entered only if KWLTB = 2 is specified in 
this namelist. NTABB must be at least 3 but no greater than 25., 

No default value is specified for NTABB. 

XWLB A one-dimensional real variable array (dimensioned at 25) consisting 

of NTABB elements. The XWLB array must be entered only if KWLTB = 2 
is specified in this namelist. Each element of XWLB specifies the 
axial (x) position, in either feet or meters, at which the wall tem- 
perature or its normal derivative is to be specified for the cowl. 
Each element of XWLB must be nonnegative, monotonically increasing, 
and within bounds of the specified geometry. No default values are 
specified for the elements of .XWLB. 

TWLB A one-dirnensional real variable array (dimensioned consisting 

of NTABB elements. The TWLB array must be entered only if KTYPE 
and KWLTB » 2 are specified in this namelist. Each element of TWLB 
specifies the wall temperature, in either R or K, at the axial posi- 
tion specified by the corresponding element of XWLB.. Each element of 
TWLB must be positive. NO default values are specified for the ele- 
ments of TWLB. 

DTDYWB a one-dimensional real variable array (dimensioned at 25) consistin 
of NTABB elements. The DTDYWB array must be entered only if KTYPE 
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2 and KWLTB = 2 are specified in this namelist. Each element of 
DTDYWB specifies the cowl wall temperature, normal derivative, in 
either (R/ft) or (K/m), at the axial position specified by the cor- 
responding element of XWLB. No default values are specified for 
the elements of DTDYWB. 


10. NAMELIST LISTS 

The parameters entered in namelist LISTS specify the wall mass bleed distri- 
bution. If the case being considered requires no mass bleed, then the parameters 
entered in this namelist do not have to be specified. For those cases which re- 
quire mass removal at the wall, the bleed distribution is specified by tabular 
input. The tabular input requires the specification of the axial location and 
extent of the bleed zones and specification of the bleed mass flux [(mass flow 
rate)/area] within each zone (see Figure 7). Options exist to specify a zero 
bleed rate for either the forebody/centerbody or cowl. 

The bleed distribution for the forebody/centerbody is specified by entering 
the following parameters. 

KDFA An integer variable denoting whether or not mass transfer occurs at 

the forebody/centerbody wall. If KDFA. = 0, a case with no mass trans- 
fer (bleed) is specified. If KDFA = 1, a case with mass transfer is 
specified. A default value of 0 is specified for KDFA. 

If KDFA = 1 is specified, the following four parameters must be entered. 

NRA A positive integer variable denoting the number of bleed zones on the 

forebody/centerbody (see Figure. 7). NRA must be at least 1 but no 
greater than 25. No default value is specified for NRA. 

X$A A one-dimensional real variable array (dimensioned at 25) consisting 

of NRA elements. Each element of XSA denotes the axial (x) position, 
in either feet or meters, of the beginning of a mass bleed zone on 
the forebody/centerbody (see Figure 7). The elements of XSA must be 
positive, monotonically increasing, and within bounds of the specified 
geometry. No default values are specified for the elements of XSA. 

XEA A. one-dimensional real variable array (dimensioned at 25) consisting 

of NRA elements. Each element of XEA denotes the axial (x) position, 
in either feet or meters, of the end of a mass bleed zone on the fore- 
body/centerbody (see Figure 7). The elements of XEA must be positive, 
monotonically increasing, and within bounds of the specified geometry. 
Koneover, for the. I^ element, XEA(I) ^XSA(I). No default values 
are specified for the elements Of XEA. 

R0VA A one-dimensional real variable array (dimensioned at 25) consisting 

of NRA elements. Each element of 2 R 0 VA denotes the mass flux, (density 
X velocity), in units of (slugAft • sec)) of (kg/(m^ • Sec)), for 
the bleed zone defined by the corresponding elements of XSA and XEA. 
For mass bleed, the elements of R 0 VA must be specified as positive. 
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FIGURE 7 BLEED DISTRIBUTION SPECIFICATION 


No default values are specified for the elements of R0VA. 

If Only the forebody flow field is to be Calculated [KCALL(l) = 1, KCALL(2) = 
KCALL(3) 0 specified in namelist LISTl], then no other parameters need be speci- 
fied in this namelist. If the internal flow is to be calculated [KCALL(2) = .1 or 
KCALL(3) = 1 specified in namelist LISTl] and mass transfer occurs at the inlet 
walls, then the following parameters must be entered to specify the cowl bleed 
distribution. 

KDFB An integer variable denoting whether or not mass transfer occurs at 

the cowl wall. If KDFB = 0, a case with no mass transfer is speci- 
fied. If KDFB = 1, a case with mass transfer (bided) is specified. 

A default value of 0 is specified for.KDFB. 

If KDFB =1 is specified, the following four parameters must be entered. 

MRB A positive integer variable denoting the number of bleed zones on the 

cowl (see Figure 7). NRB must be at least ! but no greater than 25.. 

No default value is specified for NRB. 

XSB A one-dimensional real variable array (dimensioned at 25) consisting 

of NRB elements. Each element of XSB denotes the axial (x) position, 
in either feet or meters, of the beginning of a mass bleed zone on 
the cowl (see Figure 7). The elements of XSB must be positive, mono- 
tonically increasing, and within bounds of the specified geometry. 

No default values are specified for the elements of XSB. 

XEB A one-dimensional real variable array (dimensioned at 25) consisting 

of NRB elements. Each element of XEB denotes the axial (x) position, 
in either feet or meters, of the end of a mass bleed zone on the 
cowl (see Figure 7). The elements of XEB must be positive, monotoni- 
cally increasing, and within bounds of the specified geometry. More- 
over, for the Ith element, XEB(I) > XSB(I).. No default values are 
specified for the elements of XEB... 

R0VB A one-dimenSional real variable array (dimensioned at 25) consisting 

of NRB elements. Each element of R0VB denotes the mass flux (density 
X velocity),, in units Of (slug/(ft • sec)) or (kg/(m • sec)), for 
the bleed zone defined by the corresponding elements of XSB and XEB. 
For mass bleed, the elements of R0VB must be specified as positive. 

No default values are specified for th6 e^^^ 


11. NAMELIST LIST9 

The parameters entered in namelist LISTS specify the turbulence model and 
the transition model used in the boundary layer computation. The parameters 
entered in this namelist need Only be Specified if the boundary layer computa- 
tion is to be performpd (KBLAY = 1 specified in namelist LIST!). 

The three-dimensional turbulence model incorporated into the program is 
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based on a two-layer formulation with [see. Reference (1)] 


®x " Si 

(0 ly 1 yy) 

(13) 

1 S " ^xo 

(yj .1 y 1 i) 

(14) 


where e denotes the eddy viscosity in the streamwise momentum equation, e • 
'^®9ion eddy viscosity, is the outer region eddy viscosity^^y 
is the distance measured normal to the wall, y^ is defined by the point where 
®xi ~ 6xo» §nd 6 is the local boundary layer thickness. 

The inner layer eddy viscosity is given by 





1/2 


(15) 


where u and w denote the velocity components in the boundary layer streamwise 
U) and cross-flow (z) directions, respectively, L is the mixing length, and 
ofp IS an appropriate factor accounting for the transition from laminar to tur- 
bulent flow. The mixing length L is given by 


L = K> [1 - exp (-y/A)] (16) 

where k is the von Karman parameter taken to be a constant with a recommended 
value of 

K = 0.40 (17) 


and A is defined by 


A = A'^-f (18) 


In equation (18), f is a function of the flow gradients and local properties 
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and is presented in Reference (1). The parameter is the van Driest damping 
factor wh.ich is taken to be a constant with a recommended value of 


A* = 26.0 


(19) 


The outer region eddy viscosity is given by the velocity defect relation 


So “ ^TR“ 


' 0 . 


[ (C2.w2)V2]<^-I 


( 20 ) 


where the subscript e denotes the boundary layer edge conditions, and 4 x is 
taken to be a constant with a recommended value of 


a = 0.0168 


( 21 ) 


With the streamwise momentum equation eddy viscosity ex defined by the above 
relations, the cross-flow momentum equation eddy viscosity e^ is defined by 


e 


z 



( 22 ) 


where X is a constant. For isotropic turbulence, X takes the value 


, X = 1.0 

The turbulent eddy thermal conductivity cq is given by 

e = i 
0 Prt 


(23) 


(24) 
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where is the turbulent Prandtl number which is assumed to be a constant in 
the analysis, and e is defined by 


E 



. 43I/2 


(25) 


The three-dimensional turbulence model is defined by entering the following 
parameters. 


APLUS A positive real variable denoting the van Driest damping factor in 

equation (18). The default and recommended value of APLUS is 26.0. 

XKAPPA A positive real variable denoting the von Karman parameter in equa- 

tion (16). The default and recommended value of XKAPPA is 0.40. 

XALPHA A positive real variable denoting the constant o.in equation (20), 
The default and recommended value of XALPHA is 0.0168. 

PRT A positive real variable denoting the turbulent Prandtl number in 

equation (24). The default and recommended value of PRT is 0.90. 

TFACTR A positive real variable denoting the constant X in equation (22). 
The default and recommended..value of TFACTR is 1.0. 


In addition to the above parameters, two other input parameters are employed 
in the turbulence model of the Adams algorithm which is used to internally gener- 
ate the forebody/centerbody boundary layer initial data (KBLIDA - 1 specified in 
namelist LIST6). The user should consult Reference (5) for further discussion of 
these parameters. In general, the program is executed. by retaining the following 
two parameters at their default values. 

AKLM A positive real variable used in the Adams algorithm turbulence model 

[see Reference (5) for discussion]. The default and recommended 
value for. AKLM is 0.435. 

ALAM. A positive real variable used in the Adams algorithm turbulence model 

[see Reference (5) for discussion]. The default and recommended 
value for ALAM is 0.09. 

A number of models have been incor.porated into the computer program for de- 
termining the transition parameter <Str in equations (15) and (20). The. transition 
models include: 

Model No. 1: An instantaneous transition model with fully turbulent flow 
being specified as occuring at a particular value xj of the boundary layer 
streamwise curvilinear coordinate x_(See Figure 5). the transition parameter 
67R is calculated from 
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Str • 0.0 

(x < X.J.) 

(26) 

Str= 1.0 

(x 1 Xj) 

(27) 

I'nstantaneous transition model with fully turbulent flow 
Sr Js g?ve^ occunng at a critical Reynolds number R^^ where the Reynolds 

e u 


(28) 

”^^*^.Re» “e* Ug being the boundary layer edge. density, streamwise velocitv 

-^P^«1vely. The transit, ’on paraller'':;r!‘"’ 

fijR = 0.0 

<"e < 

(29) _ 

6tr = 1.0 

0»e i Re^l 

(30) 

Model No. 3: A progressive transition model with the 
flow occunng at xji and with the onset of fully turbulent 
The transition parameter 6jr is calculated from 

onset of transitional 
flow ocCuring at xj 2 . 

'• •. ^TR " 0-0 

(x < x^^) 

(31) 


(x > Xto) 

(32) 

6tr = ^ - ^T1 


x^2 - x.|.^ 

(Xjl ^ X < Xyg) 

(33) 
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where-x Is the local value. of the boundary layer streamwlse curvilinear coordi- 
nate. 

Model No^ 4: A progressive transition model \vith the onset of transitional 
flow occuring at the Reynolds number Reji and with the onset of fully turbulent 
flow occuring at Rgj 2 * The transitional parameter 6 xr is calculated from 


6 


TR- 


0.0 


(Re " *^0 ) 

e _e^^ 


(34) 


^TR “ 


(Rfi ^ Re ) 
e e ^2 


Re -Re 


’TR 


T1 


Re ■ Re 

®T2 ®T1 


(R._. < R. < R._ ) 


*T1 


T2 


(35) 

(36) 


where Rg is defined by equation (28). 

Model No. 5: A progressive transition model using the Dhawan and Narasimha 
formulation (8). In this model the onset of transitional flow is specified as 
occuring at xji and it is assumed that fully turbulent flow occurs at x = 2><xi. 
The transition parameter 6 jr is calculated from 


^TR 


(x < 


(37) 


6 


TR 


1 .0 - exp 


— 

r " ^ 

2-| 

-3.6097 

5^ - 1.0 



*T1 

^ • J 



(x ^ 


(38) 


The transition models are specified by entering the following one-dimensional 
and two-dimensional arrays. For each array, K * 1 and K = 2 denote the forebody/ 
centerbody and cowl boundary layers, respectively. For the two-dimensional 
arrays, the index IB denotes the IB;Ui circumferential station (see Figure 6). 

ITRANM A one-dimensional integer variable array conststinq of two elements. 

ITRANM(l) and ITRANM(2) denote the transition models to be used in 
the forebody/centerbody and cowl boundary layer computations, 
respectively. ITRANM(K) can have values of 1, 2, 3, 4, and 5 cor- 
responding to the aforementioned transition models. Default values 
of 5 are specified for both elements of ITRANM. 
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XT 


RET 


XTl 


XT2 


RETl 


RET2 


A two-dimensional real variable array dimensioned at (16, 2). 

XT (IB, K) (IB » 1 MBLAY) denotes the transition distance 

(xt) for the IBth circumferential station as described in Model 
No. 1 and must be entered only If ITRANM(K) =1. No default 
values are specified for-the elements of XT. 

A one-dimensional real variable array consisting of two elements. 

The RET array must be entered only if ITRANM(K) » 2 Is specified 
in this namelist. The elements of RET specify the Reynolds number 
[defined by equation (28)] at which the flow transitions instanta- 
neously from laminar to turbulent. RET(l) and RET(2) designate 
the transition Reynolds numbers for the forebody/ centerbody anxl 
cowl, respectively. RET(2) must be entered only If the cowl 
boundary layer is to be calculated. NO default values are speci- 
fied for the elements of RET. 

A two-dimensional real variable array dimensioned at (16, 2). 

XTl (IB, K) (IB = 1, ...» MBLAY) denotes the transition distance 
(Xj,) for the IBth circumferential station as described in Model 

Noi. 3 and 5, and must be entered only if ITRANM(K) = 3 or S. No 

default values are specified for the elements of XTl. 

A two-dimensional, real variable array dimensioned at (16, 2). 

XT2.(IB, K) (IB = 1, .... MBLAY) denotes the transition distance 
(x-r^) for the IB^ circumferential station as described in 
Model No. 3, and must be entered only If ITRANM(K) * 3. No 
default values are specified for the elements of XT2. 

A one-dimensional real variable .array consisting of two elements. The 
RETl array must be entered only If ITRANM(K) 4 is specified in this 
namelist. The elements of RETl specify the Reynolds number [defined 
by equation (28)] at which the onset of transitional flow occurs. 
RETl(l) and RET1(2) specify the transitional flow onset Reynolds 
number for the forebody/centerbody and cowl, respectively. RET1(2) 
must be specified only if the cowl boundary layer is to be calculated. 
No default values are specified for the elements of RETl. 

A one-dimensional real variable array consisting of two elements. The 
RET2 array must be entered only if ITRANM(K) =4 is specified in this 
namelist. The elements of RET2 specify the Reynolds number [defined 
by equation (28)] at which the onset of fully turbulent flow occurs. 
RET2(1) and RET2(2) specify the fully turbulent flow onset Reynolds 
number for the forebody/centerbody and cowU respectively. RET2(2) 
must be Specified only if the cowl boundary layer is to be calculated. 
No default values are specified for the elements of RET2. 
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12. NAMELIST LIST10 


The parameters entered iii namelist LISTIO specify the various convergence 
tolerances and iteration limits used in the numerical integration. All para- 
meters In this namelist have specified default values. In general, the program 
is. execut£d_w1thout changing tii e values of any of the parameters in this name- 
list. 

SAFEIN A positive real variable denoting the ratio of the axial marching 
step taken to the axial marching step allowed by the Courant- 
Friedrichs-Lewy (CFL) stability criterion. This variable is used 
to determine, the axial position of both the first solution plane in 
the fore body flow field integration [KCALL(l) =l.]and the first 
solution plane in the internal flow field integration in which 
shock waves are not discretely fitted [KCALL(3) = 1]. Ensuing solu- 
tion planes for these integration options have their axial, locations 
adjusted in accord with an internally computed value of SAFEIN. 

For the internal flow field integration option in which shock waves 
are discretely fitted [KCALL(2) = 1], the axial position of each 
solution plane (except in the vicinity of a shock wave reflection) 
is controlled by the input value of SAFEIN. The specified value of 
SAFEIN must be positive, and must be less than 1.0 to satisfy the 
CFL stability criterion. The default and recommended value of 
SAFEIN is 0.975. 

CRIT A one-dimensional real variable array consisting of 18 elements. 

Each element of CRIT specifies a convergence tolerance or other 
parameter. The elements of CRIT have the following definitions 
and default values. 

CRIT(l) A positive real variable denoting the tolerance, in either feet or 
meters, used to determine if a user supplied initial-value plane 
data point is Sufficiently close to a plane of Symmetry when the 
data point is supposed to lie - the plane of symmetry. CRIT(l) is. 
also used to determine if a user supplied initial-value plane data 
point is sufficiently close to the solid boundary when that point 
is supposed to lie on the solid boundary. A default valueof 0.1 
^ ft (or 0.1m) is specified for CRIT(l). 

I CRIT(2) A positive real variable denoting the, relative tolerance used in 

I testing for a loss of significance in IBM library subroutine GELG 

1 (GELG is used to solve a system of Simultaneous linear equations). . 

' A default value of 10“’^ is specified for CRIT(2). 


1 
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CRIT(3) 


CRIT(4) 


CRIT(5) 


CRIT(6) 


CRIT(7) 


CRIT(8) 


CRIT(9) 


A positive real variable denoting the relative tolerance used in 
testing for convergence in the internal generation of the initial- 
value plane flow property field. A default value of lO"*^ is 
specified for CRIT(3). 

A positive real variable denoting the relative tolerance used in 
testing for convergence of all three coordinates in the iterative 
scheme employed in computing a streamline-surface intersection or 
a bicharacteristic-surface intersection. A default value of 10" 
is specified for CRIT(4). 

A positive real variable denoting the relative tolerance used in 
testing for the convergence of the five flow properties u, v, w, 

P, and p in subroutine S0LVE. A default. value of 10"^ is specified 
for CRIT(5). 

A positive real variable denoting the relative tolerance used in 
testing for the convergence of the static pressure in subroutine 
SH0CK. Convergence is attained in the local iteration loop if 


1P(2)-P*(2)|/P(2) < CRIT(6) 


where P(2) is the solution point pressure obtained from the local 
Hugoniot relations, and P*(2) is the pressure obtained from the 
wave surface compatibility relation. A default value of 10"'^ is 
specified for CRIT(6). 

A positive real variable denoting the mass flow rate ratio at 
which the maximum number of radial stations, allowed in the fore- 
body flow field computation is changed from JLIMIT(I) to JLIMIT(2). 
The mass flow rate ratio is the mass flow rate at a given fore- 
body flow field solution plane divided by the estimated mass flow 
rate at the axial station corresponding to XEND(l). A default 
value o-f 0.5 is specified for CRIT{7). 

A positive real variable used as a multiplier of the mass flow 
ratio which is employed in determining whether or not point 
addition is to be performed on a solution plane in the forebody 
flow fie.ld integration. The mass flow rate ratio is the mass flow 
rate at the solution plane just computed divided by that at the 
last solution plane where point addition or deletion was performed. 
A default value of 1.0 is specified for CRIT(8). 

A positive real variable denoting the relative tolerance used in 
routine LINK31 for determining when the angle a calculated in the 
global correction for the bow shock wave points has converged.. A 
default value of 10“4 is specified for CRIT(9). 
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CRIT(IO) 


CRIT(ll) 


CRIT(12) 


CRIT(13) 

CRIT(14) 


A positive real variable used in routine LINK31 for determining if 
a sufficient numbdr of shock wave solution points have converged in 
global correction. Convergence is attained when 


M/IST0P > CRIT(IO) 


where M is the number of shock wave solution points which have 
converged in global correction, and IST0P is the number of cir- 
cumferential stations in the computed sector. A default value 
of 0.8 is specified for CRIT(IO). 

A positive real variable denoting the relative tolerance used in 
subroutine BSH0CK for determining if the velocity component down- 
stream of the shock wave and normal to the surface of the solid 
boundary has converged to a specified value, A default value of 
10"^ is specified for CRIT(ll). 

A positive real variable used in routine LINK32 for determining if 
another solution plane is to be inserted between the last solution 
plane and the intersection of the incident internal shock wave with 
the solid boundary. Another solution plane is inserted if 


Ax/AX(.pL > CRIT(12) 


where ax is the axial (x) distance between the last computed plane 
and the nearest point on the space curve defined by the intersec- 
tion of the incid-int interanl shock wave with solid boundary, and 
Axrci is the axial step allowed by the Courant-Friedrichs-Lewy (CFL) 
stability criterion. A default value of 0.2 is specified for CRIT(12). 

Not presently employed. 

A positive real variable used in subroutine PENTRE for determining 
if a streamline-shock wave intersection point is sufficiently close 
to the current solution plane, so that an interior point unit process 
on the downstream side of the shock wave is not performed. Instead, 
a streamline projection onto the solution plane and subsequent flow 
property interpolation in this plane is performed. The application 
of the interior point unit process is not performed if 


(x^ - ^ GRIT(14) 
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CRIT(15) 


CRIT(16) 


CR1T(17) 


CRIT(18) 


where x< is the axial position of the solution plane, x^nt 

axial location of the streamline.-shock wave intersection point,, and 

AxcFL "is the axial marching step allowed by the 

Lewy stability criterion. A default value of 0.4 is specified fo 

CRIT(14). 

A oositive real variable used in. subroutine STRSHK for. determining 
if convergence has been obtained in calculating. the intersection 
point of a body streamline with the space curve defined by the 
intersection of the incident internal shock wave with a solid 
boundary. Convergence is attained when 


:^1+1 ■“ V 


CRIT(15) 


where 6i is the polar angle of the 

iteration, and ©i+i is the polar angle on the 

A default value of 10“4 radians is specified for CRIUlb)- 

A positive real variable used in subroutine INTSCT for determining 
if convergence has been obtained in calculating the intersection 
point of a bicharacteristic with either a solid boundary or a shock 
wave, or the intersection point of a streamline with a shock wave. 
Convergence is attained when 


iRl - Rgl -- CR^T(16) 


where Ri is the radius of the intersection point obtained by inte- 
grating the equation for a streamline or bicharacteristic, and Rg 
is the intersection point radius obtained from the shock vvaye or 
boundary surface formulations. A default value of 10" ft tor 
lir‘*m) is specified for CRII(16). 


A positive real variable denoting the relative tolerance used in 
subroutines ABLSLN and GBLSLN to test for convergence of the wall 
boundary streamwise velocity normal derivative in the boundary layer 
implicit finite difference algorithm. A default value of 10" is 
specified for CRIT(17). 

A oositive real variable denoting the tolerance used in subroutine 
SBLINT to test for convergence of the continuity, streamwiSe momen- 
tum, and cro.«vS-flow momentum equations in the shock. wave-boundap 
layer integral- interaction analysis. A default value of 10 is 
specified for CRiT(18). 


ITEND A one-dimens1oi1al integer variable array consisting of 8 elements. 

Each element of ITEND specifies a limit to the number of iterations 
permissible in a given iteration loop. The elements of ITEND have 
the following definitions and default values. 

ITEND(l) A positive integer variable denoting the maximum number of inner 

interations permissible in determiniJig the intersection coordinates 
of either a streamline with a surface, or a bi characteristic with a 

surface. ITEND(l) is used in conjunction with C^IT(4). A default 

value of 10 is specified for ITEND(l). 

ITEND{2) A positive integer variable denoting the maximum number of outer 
iterations permissible in obtaining the five flow properties u, v, 
w, P, and p in all unit processes except the shock wave-solid boundary 
point unit process. A default value of 10 is specified.for ITEMD(2). 

ITEND(3) Not presently employed. 

ITEND(4) A positive integer variable denoting the maximum number of iterations 
permissible in the relaxation of the velocity component normal, to the 
solid boundary and downstream of the reflected (cowl lip) shock wave 
in the shock wave-solid boundary point unit. process (subroutine 
BSH0CK). A default value of 20 is specified for ITEND(4). 

ITEND(5) Not presently employed. 

ITEND(6) A positive integer variable denoting the maximum number of permissible 
subiterations in determining the intersection point of a line segment 
with a given three-dimensional surface (subroutine INTSCT). A de- 
faull value of 10 is specified for ITEND(6).,. 

ITEND(7) A positive integer variable denoting the maximum number of itera- 
tions permissible in obtaining convergence for the wall boundary 
normal derivative of the streamwise velocity in the boundary layer 
implicit finite difference algorithm. ITEND(7) is used in conjunc- 
tion with CRIT(17). A default „Vja.lue oflOis specified fOr ITEND(7). 

ITEND(8) A positive integer variable denoting the maximum number of iterations 
permissible in obtaining convergence for the continuity, streamwise 
momentum, and cross-flow momentum equations in the integral shock 
wave-boundary layer interaction analysis. ITEND(.8) is used in 
conjunction. with CRIT(18). A default value of 15 is specified for 
ITEND(8). 


13. NAMELIST LISTll 

The parameters entered in namelist LISTll specify if debug output is to 
be printed. 


KDUMR-. A one-dimensiOnal integer variable array consisting of 9 elements. 
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Each element of KDUMP specifies whether or not a particular -compu- 
tational process is. to have debug output printed. Specifying 
KDUMP(I) = 1 (1=1 to 9) activates the debug, output option for the 
corresponding process. Specifying KDUMP{I) = 0 causes no debug 
output to be printed for the Corresponding process. The elements 
of KDUMP activate the debug Output option for the following processes 


and have 

the_following default values. 


KDUMP( I ) 

Activates Debug Output for 

Default Value 

KDUMP(l) 

supersonic flow interior 
point scheme 

0 

KDUMP(2) 

supersonic flow solid body 
point scheme 

0 

KDUMP(3) 

supersonic flow field-shock 
wave point scheme 

0 

KDUMP(4) 

supersonic flow solid body- 
shock wave point scheme 

0 

KDUMP(5) 

boundary layer initial data 
generation schemes 

0 

KDUMP(6) 

interpolation scheme for 
boundary layer edge conditions 

0 

KDUMP(7) 

plane of symmetry finite dif- 
ference -boundary layer scheme 

0 

KDUMP(8) 

three-dimensional finite dif- 
ference boundary layer scheme 

0 

KDUMP(9) 

shock wave-boundary layer 
interaction region scheme 

0 


ISTART A positive integer variable denoting the solution plane number at 
which debug output is to be initiated. A default value of 1 is 
specified for ISTART. 


14. .FORMATTED READ OF THE INITIAL-VALUE FLOW PROPERTY FIELDS 

The user supplied supersonic flow initial data are entered by a formatted 
read of file ITAPl after all eleven namelists have been input. To enter the 
supersonic, flow initial data_by tabular input, KIVS = 0 must be specified in 
namelist LIST2. The. default value for ITAPl is 5 (the Input file). The super- 
sonic flow initial-value plane may be read from TAPEIO by specifying ITAPl. = 

10 in namelist LIST2. TAPEIO is linked to the dummy file IVSl in the PR06RAM 
card. 


The index limits for the supersonic flow initial-value plane, IST0P and 
JMAXI, are entered in namelist LIST3. The initial-value plane point networks 
for the four flow tymmetry options are illustrated in Figure 3.. 
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The supersonic flow Initial data are entered by the formatted read state* 

ment 

READ (ITAPl.II.l) ((Y(IiJ),Z(I,J.),U(I,J),V(I,J),W{I,O),P(I,J),R0(I,j), 

J-1,JMAXI),I=1,IST0P) 

with the format (4E20.13/3E20.13) . The parameters In the formatted read state- 
ment have the following definitions (see Figure 3). 

I - An integer denoting the circumferential Index of the data point. 

J An integer denoting the radial -Index of the data point. 

Y A two-dimensional real variable array, each element of which denotes 

the y-positlon, in either feet or meters, of point The Y 

array Is dimensioned at (30,15), No default values are specified 
for the elements of Y. 

Z A two-dimensional real variable array, each element of v^hich de- 

notes the 2-pos1t1on, In either feet or meters, of point (I,J). 

The Z array Is dimensioned at (30,15). No default values are speci- 
fied for the elements of Z. 

U A two-dimensional real variable array, each element of which de- 

notes the x-component of velocity, in either (ft/sec) or (m/s), at 
point (1,0). The U array Is dimensioned at (30,15). No default 
values are specified for the elements of U. 

V A two-dimensional real variable array, each element of which denotes 
the y-component of velocity, in either (ft/sec) or (m/s), at point 
(1,0). The V array 1s dimensioned at (30,15). No default values 
are specified for the elements of. V. 

W A two-dimensional real variable array, each element of which denotes 

the 2-component of the velocity, in either (ft/sec) or (m/s), at 
point (1,0). The W array Is dimensioned at (30,15). No default- 
values are specified for the elements of W. 

P A two-dimensional real variable array, each element of which denotes 

the pressure, in either. (Ibf/ft2) or (N/m2), at point (1,0). The 
P array is dimensioned at (30,15). No default values are specified 
for the .elements of P. 

A two-dimensional real variable, array, each element of -which denotes 
the density, In either (slug/ft^) or (kg/m3), at point (1,0). The 
R0 array Is dimensioned at (30,15). No default values are specified 
for_the elements of R0. 

In all Cases, the supersonic flow_1n1t1al-value plane data points with 0 = 1 
must lie (to a close approximation) on the surface of the forebody/centerbody. 

The supersonic flow Initial-value plane data paints with J = JMAXI correspond to 
the downstream bow shock wave points. If only the internal flow. Is to be computed, 
the Initial-value plane is located at the cowl. Up axial station. It is Suffi- 
cient to specify the flow property field, in this case, to a point just outside 
of the Cowl lip. If the bow shock wave radius is less than that of the Cowl 
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lip, the execution is aborted. 


For the case of no planes of flow symmetry (KSYM = 0), the computed sector 
is the entire solution plane corresponding to IST0P circumferential stations 
and JMAXI radial stations [see Figure 3(a)]. For the case of one plane of flow 
symmetry (KSYM =1), the computed sector is the half-plane bounded by the y-axis 
and containing the +z-axis.[see Figure 3(b)]. In this case, the data points 
with I = 1 must lie on the +y-axis, and the data points with I = IST0P must lie 

on the -y-axis. For the case of two planes of flow symmetry (KSYM = 2). the com- 

puted sector is the quadrant bounded by the +y-axis and the +z-axis [see Figure 
3(c)]. In this case, the data points with I = 1 must lie on the +y-axis, and _ 

the data points with I = IST0P must lie on the +z-axis. For the axisymmetric 

flow case (KSYM = 3), the computed sector is limited to the single circumferen- 
tial station lying on the +y-ax1s [see Figure 3(d)]. In this case, the data 
points with I = 1 must lie on the +y-axis. 

If the forebody supersonic flow field is not being calculated [KCALL(l) = 0], 
or if it is being computed and the bow shock wave is conical (KC0N =1), then the 
following shock wave angle input parameters do not have to be entered. If, how- 
ever, the forebody flow field is to be calculated [KCALL(l) * 1] and the bow shock 
wave is not conical (KC0N - 0), then the angle subtended by the bow shock wave 
and the x-axis in the meridional plane defined by the shock wave point for each 
shock wave point In the computed sector must be entered by the foramtted read 
statement. 


READ (ITAP1,II2) (BETA( I) ,1=1 ,IST0P) 


with the format (E20.13). The parameters in the formatted read statement have 
the following definitions (see Figure 3). 

I An integer denoting the circumferential index of the initial-value 

plane shock wave point. 

BETA A one-dimensional real variable array (dimensioned at 30), each 

element of which denotes the angle, in radians, subtended by the 
bow shock wave and the x-axis in the meridional plane defined by 
the corresponding initial-value plane downstream shock wave point. 

Each element of BETA must be positive. No default values arO 
specified for the elements of BETA. 

The following parameters must-be entered only if the boundary layer compu- 
tation is to be performed -(KDLAY = 1 specified in namelist LISTl). 

The user-supplied forebody/centerbody boundary layer initial data are entered 
by a formatted read Of file ITAP2 after all eleven namelists have been input and 
after the supersonic flow initial data have been input if that data was externally 
generated. To enter the forebody/centerbOdy boundary layer initial data by tabular 
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must be specified in namelist LIST6. The default val^e for 
UAP2 1 S 5 (the input file). The forebody/centerbody boundary layer Initial 

from TAPEll by specifying ITAP2 * 11 in namelist LIST6. TAPEll 
IS linked to the dummy file IVS2 in the PR06RAM card. 

MRi forebody/centerbody boundary layer initial data. 

3M Illustrated irpiSirre.'" boundary layer point networks 

matted reld''stat^ments''’‘’°'‘^ boundary layer Initial data are entered by the for- 

REAO 0™K,II3M(ZZ(IB)^A(IB,JB),UA(IB.JB),WA(IB.OB).XLAMDACIB.JB), 

with the foraats (4E20.13/E20.13) and (4E20.13/4E20.13), respectively. The para- 
Mters in the formatted read statements have the following definitions (see Figure 


IB 

JB 

11 


FA 


UA 


WA 


XLAfiDA 


UEA 


An integer denoting the circumferential index of the data point. 

An integer denoting the radial (or boundary layer normal) index of 
the solution point. 

A one-dimensional real variable array, each element of which denotes 
the circumferential angle, in radians, of the IB th circumferential 
station. The l 1 array is dimensioned at (16). No default values 
are specified for the elements of 22. 

A two-dimensional real variable array, each element of which denotes, 
the distance, in either feet or meters, between point (IB,JB) and 
^ +^rTD*'1n\^ measured along the surface normal passing through 
point IIB.JB). The FA array is dimensioned at (16,20). No default 
values are specified for the elements of FA. 

variable array, each element of- which denotes 
the boundary layer streamwise(x) component of velocity, in either 
(ft/sec) or (m/s), at point (IB,JB). The. UA array is dimensioned 
at \ib,20).. No default values, are specified for the elements of UA. 

A two-dimensional real variable array, each element of which denotes 
the boundary layer cross flow (z) component of velocity, in either 
(ft/sec) or (m/s), at point (IB,JB). The WA is dimensioned at (16,20) 
No default values are specified for the elements of WA. 

A two-dimensional real- variable array, each element of which denotes . 

either (slug/ft3) or (kg/m^), at point (IB,JB). The 
XLAMDA array is dimensioned at (16,20). NO default values are speci- 
fied for the elements of XLAMDA. 

A One-dimensional real variable array, each element of which denotes 
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the boundary layer edge streamwise (x) component of velocity, in 
either, (ft/sec) or (m/sec), at the circumferential station. 

Thd UEA array is dimensioned at (16). No default values are speci- 
fied for the elements of UEA. 

WEA A one-dimensional real variable array, each element of which denotes 

the boundary layer edge cross flow (z) component Of velocity, in either 
(ft/seC) Or (m/sec), at the IBtl]^ circumferential station. The WEA 
array is dimensioned at (16). No default values are specified for. 
the elements of WEA. 

PEA A one-dimensional real variable array, each element Of which denotes 

the boundary layer edge pressure, in either (Ibf/ft^) or (N/m^), at 
the IB^ circumferential station. The PEA array is dimensioned. at 
(16). No default values are specified for the elements of PEA. 

R0EA A one-dimensional real variable array, each element of which denotes 

the boundary layer edge density, in either (slug/ft3) or (kg/m^), 
at the IB^ circumferential station. The R0EA array is dimensioned 
at (16). No default values are specified for the elements of R0EA. 

DUDXA A one-dimensional real variable array, each element of which denotes 

the boundary layer edge derivative of the streamwise velocity compo- 
nent (ue) with respect to the streamwise curvilinear coordinate (x) 
at the IBUi circumferential station. The units for DUDXA are (sec-i). 
The DUDXA array is dimensioned at (16). No default values are speci- 
fied for the elements of DUDXA. 

DWDXA A one-dimensional real variable array, each element of which denotes 

the boundary layer edge derivative of the cross flow velocity compo- 
nent (wg) with respect to the streamwise curvilinear coordinate (x) 
at the IB-^ circumferential station. The units for DWDXA are (sec~i). 
The DWDXA array is dimensioned at (16). No default values are speci- 
fied for the elements of DWDXA. 

DPDXA A one-dimensional real variable array, each element of which denotes 

the boundary layer edge derivative of the static pressure (Pg) with . 
respect to the streamwise curvilinear coordinate (x) at the IB^ 
circumferential station. The units for DPDXA are either (lbf/ft3) or 
(N/m3). The DPDXA array is dimensioned at (16). No default values 
are specified for the elements of DPDXA. 

DRDXA A one-dimensional real variable array, each element of which denotes 

the boundary layer edge derivative of the density (pe) with respect 
to the streamwise curvilinear coordiante (x) at the IB^ circumfer- 
ential station. The units for DRDXA are either (slug/ft^) or (kg/m^). 
The DRDXA array is dimensioned at (16). No default values are speci-- 
fied for the elements of DRDXA. 

The user-supplied cowl boundary layer initial data are entered by a formatted 
read of file ITAP3. after all eleven namelists have been input and after the super- 
sonic flow and forebody/centerbody boundary layer flow initial data have been input 
if they are externally generated. To enter the cowl boundary layer initial data 
by tabular input, KBLIDB = 0 must be Specified in namelist LIST6. The default 
value for 1TAP3 iS 5 (the input file). The Cowl boundary layer initial data may 
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b6 read from TAPE12 by specifying ITAP3 = 12 in namelist LISTS. TAPE12 is 
linked to the dummy file IVS3 in the PROGRAM card. 

The index limits for the cOwl boundary Uyer initial data, MBLAY and NINB, 
are entered in namelist LISTS. The boundary layer point networks are illustrated 
in Figure 6. 

The cowl boundary layer initial data are entered by the formatted read state 

ments 


READ {ITAP3,II3) ((ZZ(IB) ,FA(IB,JB) ,UA(IB,JB) ,WA(IB,JB) ,XLAMDA(IB,OB) , 
JB=1, NINB), IB=1, MBLAY) 

READ (ITAP3,II4) (UEA ( IB) ,WEA(I^,PEA(IB) ,R0EA(IB) ,DUDXA( IB) ,DWDXA( IB) , 
DPDXA(IB),DRDXA(IB),IB=1,MBLAY) 

with the same formats and parameter definitions as presented before for the 
forebody/centerbody boundary layer initial data formatted read. 

In all cases, the boundary layer initial data points with OB = 1 must lie 
(to a close approximation) on the appropriate body surface. For the case of 
one plane of flow symmetry (KSYM = 1), the initial data points with IB = 1 must 
lie on the -y-axis, and the points with IB = MBLAY must lie on the +y-axis. 
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SECTION V 


OUTPUT INTERPRETATION 

1 . INTRODUCTION 

The initial portion of the computer output comprises preliminary infor- 
mation. This preliminary output consists of information identifying the 
problem being considered, the specified computation options, the flow symmetry 
option, the thermodynamic model and the molecular transport properties, the 
vehicle orientation and the free-stream conditions, certain index parameters, 
the contours of the centerbody and the cowl, and the convergence tolerances 
and the iteration limits. The initial -value plane is then printed. Alter- 
natively, if a program restart is specified, the last solution plane written 
on the restart file is printed. Each solution plane is then printed in a 
format similar to the initial-value plane printout. The supersonic core flow 
solution is printed first, followed by the forebody/centerbody and cowl 
boundary layer solutions. Additionally, the redistributed data plane at the 
cowl lip axial station is printed if the internal flow integration option is 
specified. Moreover, for the internal flow field computation, the solution 
points. are printed which lie along the space curves defined by the intersection 
of the internal shock wave with the solid boundaries. 


2. SUPERSONIC CORE FLOW SOLUTION OUTPUT . I 

The supersonic core flow solution output parameters listed on the computer 
printout are defined below. 1 

I circumferential index of the Solution point ; 

-♦ 

J radial index of the solution point j 

X axial position of the solution plane or the solution point, (ft) or (m) ■ 

4 

Y y-position, (ft) or (m) ‘ 

Z z-position, (ft) or (m) | 

M Mach number i 

Q velocity magnitude, (ft/$ec) or (m/sec) 

P pressure, (Ibf/ft^) or (N/m^) 



R0 density, (slug/ft^) or (kg/iii^) 

T absolute temperature, (R) or (K) 

U x-component of velocity, (ft/sec) or (m/sec) 

V y-component of velocity, (ft/sec) or (m/sec) 

W z- component of velocity, (ft/Sec) or (m/sec) 

p 2 

PT Stagnation pressure, (Ibf/ft ) or (N/m ) 

TT stagnation temperature, (R) or (K) 

IT6 number of global corrector applications 
ITL number of local iterations 


Streamline solution points have both I and J indices which are numbers. 

An upstream shock wave solution point is denoted by a numerical I index, and 
the 0 index is U. A downstream shock wave solution point is denoted by a 
numerical I index and the 0 index is D. 

for the external flow field about the forebody, the body streamline 
solution points are denoted by 0 - 1. The outer bound to the computational 
flow regime is defined by the locus of downstream shock wave solution points, 
points with J = D. Since periodic point addition and deletion are performed 
in the external flow field integration, continuous streamlines throughout 
the computational flow regime are .not available. Inserted solution points 
are noted by IT6 = ITL =0. 

For the continuous internal flow field integration option, the body 
streamline points on the surface of the centerbody are denoted by J =-l, and 
the body streamline points on the surface of the cowl are denoted J = OINLET. 
The solution is found on the continuous streamlines which pass through the 
redistributed points on the solution. plane at the cowl lip axial station. 

For the internal flow field integration option in which shock waves are - 
discretely fitted, the body streamline points on the surface of the center- 
body are denoted by 0=1. The body streamline points on the surface of the 
cowl are denoted by 0 = (OINLET - 2). The shock wave solution points float 
in the storage arrays as the internal shock wave travels between. the center- 
body and the cowl on successive solution planes. The streamline points 
between the upstream side solid boundary (either centerbody or cowl) and the 
upstream Shock wave points on a given solution plane lie in the upstream flow 
field sector On that solution plane. In a like manner, streamline points 
which lie be.tween the other solid boundary and the downstream shock wave 
points on a given solution plane lie in the downstream flow field Sector on 
that solution plane. A reversal of the upstfeam and downstream sectors occurs 
at an internal shock wave-solid boundary intersection. It should be noted 
that continuous streamlines are followed in the internal flow field integra- 
tion..- 
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The -Intersection of the. incident internal shock wave with a solid bound- 
ary at a shock wave-sol id- boundary intersection defines a sp'-ce curve. The 
solution is found on both the upstream and downstream sides of both the in- 
cident and reflected shock waves at points on -this space curve. 

At the end of a solution plane printout, the Courant number and the 
x-step regulation parameters are printed. The Courant number is the ratio 
of the axial step taken to the axial step allowed by the Courant-Friedrichs- 
Lewy stability criterion (based On immediate neighbors in the interpolation 
fit. point stencils). The Courant number listed is that used in obtaining 
the solution plane-that was just printed. Likewise, the x-step regulation 
parameters refer to the solution plane that was just printed. 


3. BOUNDARY LAYER FLOW SOLUTION OUTPUT 

prmtJSt »5ted on the compute 


I . cicrumferential index of the solution point 
0 radial index of the solution point 

X axial position of the corresponding superSon'c core flow solution 
plane for...the current boundary layer .solution surface, (ft) or (m) 

XC streamwise curvilinear coordinate (x) of the boundary layer solution 
surface (for the forebody/centerbody , XC=0.0 corresponds to the fore- 
body tip; for the cowl, XC=0.0 corresponds to the Cowl lip), (ft) or 
(m) 

Y distance measured normal to the. wall, (ft) or (m) 

Z polar angle (see . Figure 6), (2=0.0° corresponds to the windward 
meridian; 2=180.0° corresponds to the leeward meridian), (degrees) 

M Mach number 

Q velocity magnitude, (ft/sec) or ...(m/sec) 

P pressure, (Ibf/ft^) or (N/m^) 

R0_ density, (Slug/ft^) or (kg/m^). 

T absolute temperature, (R) or (K) 
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U streamwise (x) coiiiponGnt of velocity, (ft/sec) or (m/sec) 

V normal (y) component of velocity, (ft/$ec) or (m/sec^ 

W cross*flow (z) component of velocity, (ft/s6c) ojL(m/$ec) 

PT stagnation pressure, (Ibf/ft^) Or (N/ni^) 

TT stagnation temperatur. , (R) or (K) 

UE streamwise (x) component of velocity at the boundary layer edge, 

(ft/sec)„ or (m/sec) „ 

WE cross-flow ( 2 ) component of velocity at the boundary layer edge, 

(ft/sec) or (m/sec) . 

PE pressure at. the boundary layer edge, (Ibf/ft^) or (N/m^) 

R0E density at the boundary layer edge, (slug/ft^) or (kg/m^) 

TE absolute temperature at the boundary layer edge, (R) or (K) 

HTE total enthalpy at the boundary layer, edge, (fAsec^) or (m^/sec^) 

DUEOX first partial derivative of the boundary layer edge streamwise 

velocity component (UE) with respect to the streamwise curvilinear 
coordinate (x), (sec'^) 

DWEDX first partial derivative of the boundary layer edge cross-flow. 

velocity component (WE) with respect to the streamwise curvilinear 
coordinate (x), (sec“^) 

DPEDX first partial derivative of the boundary layer edge pressure (PE) with 

respect to the streamwise curvilinear coordinate (x), (Ibf/ft^) or (N/m3) 

DR15EDX first partial derivative of the boundary layer edge derisity (Rf!E) 

with respect to the streamwise curvilinear coordinate (x), (slug/ft^) 
or (kg/m'^) 

DUEDZ first partial derivative of the boundary layer e.dge streamwise velocity 
component (UE) with respect to polar angle-, (ft/sec/ radian) or 
(m/sec/ radian) 
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DWEDZ first partial derivative of the boundary layer edge cross^flow velocity 
component (WE) with respect to polar angle, (ft/sec/radian) or 
(m/sec/ radian). 

DPEDZ first partial derivative of the boundary layer edge pressure (PE) 
with respect to polar angle, (Ibf/f A radian) or (N/m^/radian) , 

DR0EDZ first partial derivative of the boundary layer edge density (R0E) 

with respect to polar angle, (slug/ttVi^adian) or (kg/m /radian) 

DTX s-trcamwise Jjoundary layer displacement thickness (fi^j^)* or (m) 

DTZ cross-fJow boundary layer displacement thickness (6^^), (ft) or (m) 

MIX streamwise boundary layer momentum thickness or (m) 

MTZ cross-flow boundary layer momentum thickness (5,^2^ * 

TWX Streamwise wall shear stress component (Ibf/ft ) or (N/m ) 

2 2 

TWZ cross-flow wall shear stress component (Ibf/ft ) or (N/m ) 

TWT total wall shear stress (\^)» (Ibf/ft^) or (N/m^) 

R0VB mass bleed flux (R0VB<O denotes mass bleed), (slug/(ft *sec)) or 

(kg/(m^*sec)) 


The respective displacement and momentum thicknesses are defined by: 



(39) 


(1 - -^ ) 0 / 
P^We 


'tz 


mx 


(1 - ^ ) dy 


e e 


(40) 


O 


(41) 


(42) 


1112 


0 - :^) dy 
w. 


Pe^e 


components in the boundary layer stream* 


^md so5u?i”n“o?nJs“u“fS"-'’SA"5'. “ ' correspond to the well surfece. 

^nT,; JStU^?h?cl"e%e^" 

dimrtream“ofTsh2!rwl;t“^ 1 tedllteir'' 

(5 not computed and is 2Htp2t as ofo ~ component of velocity 
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SECTION VI 
SAMPLE CASES 


1. INTRODUCTION 

Four sample cas6s are presented in this section to illustrate the applica- 
tion of the computer program for calculating the flow field in supersonic mixed- 
compression aircraft inlets. For each of the four sample cases, a discussion 
of the problem is given, the required input data are presented, and selected 
portions of the computer output are listed. The input parameter discussions 
follow the order in which the input parameters are presented in Section IV. 

Sample Case No. 1 considers the computation of the supersonic external 
flow about the forebody of a typical mixed- compression inlet at angle of attack. . 
Sample Case No. 2 is concerned with the -computation of both the supersonic ex- 
ternal flow and the boundary layer flow for the forebody geometry considered 
in Sample Case No. 1.^ Sample Case No. 3 is concerned with the computation of 
the internal supersonic core flow for the Mach 3.5 inlet documented in Reference 
(9) at angle of attack. Sample Case No. 4 considers the computation of both the 
internal supersonic core flow and the boundary layer flow for the mixed- compres- 
sion inlet considered in Sample Case No. .3. 

It should be noted that additional sample cases may be found in Reference (3). 
Although Reference (3) does not discuss the boundary layer computational proce- 
dures, the discussion concerning the computation of the supersonic core flow is 
similar to that presented herein. 


2. SAMPLE CASE NO.. 1 

This sample Case is concerned with the computation of the supersonic external 
flow field about the axi symtnetri c forebody of a typical mixed-compression, inlet 
at incidence. The supersonic initial-value plane is generated internally in the 
program using the Jones algorithm described in Section II. 

The data deck for Sample Case No. 1 is presented in Figure 8. The first 
card of the data deck is the title card. English units are used, so KUNIT 
retains its default value of 1 in namelist LISTl. Since only the forebody flow 
field is to be computed, KCALL(2)'*0 is specified in namelist LISTl, while 
KCALL(l) and KCALL(3) retain their default values of 1 and 0, respectively. 

The forebody flow integration termination point, denoted by XEND(l), is 2.0 ft, 
the default value. Since one plane of flow symmetry exists, KSYM is retained 
at its default value of 1. The molecular transport terms are not to be. included 
in the computation, and global correction is to be performed on the bow shock 
wave solution points; henCe, KVISCY and KSCL0B are retained at their default 
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values of 0 and 1, respectively. Mass transfer effects are. not to be included 
in the computation, hence KTRANS-0 is specified in namelist LIST!. For this 
sample case, the boundary layer computation will not be performed, thus KBLAY»0 
is specified. The default value of RCAVG*0.8 ft is used for estimating the mass 
flow rate downstream of the bow shock wave at the cowl lip axial station. 

KPRINT retains its default value of 1, consequently all solution points are 
printed. IPRSTP and KSTART are both kept at their default values of 0, thus 
the execution is not terminated at a specified solution plane, nor jire any 
restart file operations performed. 

All input parameters in namelist LIST2 except for the free stream Mach 
number, MFS, retain their default values. Consequen+-ly, the free-stream pressure 
PFS and the free-stream density R0FS have values of 242.2 (Ibf/ft^) and 
0.0003622 (slug/ft^), respectively. The free-stream Mach number for this case 
is 2.5.; hence MFS-2.5. PITCH and YAW retain their default values of 1.0 and 
0.0 degrees, respectively.. Thus, the specified angle of attack is 1.0 degree. 

The axial qjosition of the initial -value plane, specified by the default value 
of XI, is 1.0 ft. Since KIVS and KC0N are both 1, the initial -value plane is 
internally generated and the bow shock wave is assumed to be conical (the fore- 
body is conical).. Since KSUPER retains its default value of 2, the supersonic 
flow initial data are generated using the Jones algorithm described in Section II. 
The parameter ITAPl is not employed since the initial-value plane is generated in- 
ternally (KIVS=1). 

All input parameters in namelist LISTS retain their default values. Conse- 
quently, IST0P and JMAXI have values of 15 and 11, resepctively, so that 15 cir- 
cumferential stations are employed, and 11 radial stations on the initial-value 
plane are specified. JLIMIT(l) and JLIMIT(2) retain their default values of 
11 and 15 radial stations, respectively. 

All input parameters in namelist LIST4 retain their default values. Thus, 
the specific heat ratio and gas constant, specified by GAMMA and R, respectively, 
have values of 1.4 and 1716.16116 (ft-lbf)/(slug-R) , respectively. Since. the. 
molecular transport terms are not included in the computation (KVISCY=0), the 
input parameters VIS0, T0, B and PR are not employed. 

All input parameters in namelist LIST5 retain their default values. The 
default inlet geometry has a conical forebody/centerbody with a cone half-angle . 
of 10.0 degrees. The forebody tip is located at x=0.0 ft, and the centerbody 
geometry is specified to x=3.5 ft. Thus, NCENT=2, KDCENT(1)*'3, XCENT(1)=1.0, 
XCENT(2)-3.5, and C0NE=1O.O. The remaining input parameters in namelist 
LIST5 are not employed. 

Since the boundary layer computation is not invoked and macs transfer 
effects are not to be considered (KBLAY»0 and KTRANS-0 specifiud in namelist 
LISTl), the input parameters in namelists LIST6, LIST7, LISTS, and LIST9 do 
not have. to be entered. 

All Convergence tolerances and iteration limits retain their default values 
in namelist LISTl 0. 
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No debug output is to be printed, hence all input parameters in namelist 
LISTll retain their default values. 

Selected 4 )ortionS of the computer output for this sample case are presented 
in Figure 9. The first portion of the program output presents the job title, 
the specified computation options, the flow symmetry option, the thermodynamic 
model, the vehicle orientation and free-stream data, the type of initial-value 
plane, certain index parameters, the centerbody and the cowl contours, and the 
various convergence tolerances and iteration limits. The next portion presents 
the internally generated initial-value plane flow property field. The final 
portion of the program output presents selected supersonic flow solution planes. 
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Figure 8. Data deck for Sample Ca^e No 
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Figure 9. Continued. 
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3. SAMPLE CASE NO. 2 


This sample-case -is concerned with the computation of both the supersonic 
external flow and ..the boundary layer flow for the forebody geometry considered 
in Sample case. No. 1. Both the supersonic flow initial, data and the boundary 
layer flow initial data are generated internal ly_in the program. 

The data deck for Sample Case-No. 2 is presented in Figure 10. Again, the 
first card of the data deck is the title card. 

Namelist LISTI for Sample Case. No. 2 is. identical to that for Sample Case 
No. 1 except for the value of KBLAY. Since in the present case, the boundary 
layer is to be computed, KBLAY is retained at its default value of 1. 

Namelists LIST2,. LISTS,. LIST4, and LISTS for Sample Case No. 2 are identical 
to those specified in Sample Case. No. .1. 

All input parameters in namelist LISTS retain their default values except 
for FACT0R and CLENCH. As a consequence, KBLIDA and KTURB retain their default 
values of 1 and 0, respectively. Thereby, the forebody boundary layer initial 
data are generated internally and laminar flow is specified.. CLENCH is specified 
as 3.5 ft, the value to which the forebody geometry is specified. The Adams 
program input parameters ADY and ARATI0 are retained at their default values of 
0.010 and 1.0630, respectively. Since the boundary layer initial data aro being 
internally generated, the input. parameter XASTRT, which denotes the position of 
the initial data surface, does not have tn bo entered but rather is computed 
internally. Likewise, the input parameter ITAP2 is not specified.. MBLAY and 
NA are retained at their default values of 15 and 20, respectively. Hence, 

15 circumferential and 20 radial stations are used in the forebody boundary 
layer calculation. The boundary layer normal coordinate stretching factor is 
selected as 1.10; hence, FACT0R(1)=1 . 10 is specified in namelist LISTS. Note 
that since the cowl boundary layer computation is not invoked, the cowl boundary 
layer.input parameters do not have to be specified. 

All input parameters in namelist LIST? retain their default values except 
for TC0NS.TA. Consequently, KTYPE retains its default value of 1, which speci- 
fies that a wall temperature rather than temperature derivative boundary condi- 
tion will be used in the boundary layer computation. The input parameter, 

KWLTA retains its default of 1, hence a constant wall temperature boundary 
condition.wi.il be employed. Since KWLTA is specified, the wall temperature is 
specified by entering TC0NSTA, which has been selected as 876.7 R, the free- 
Stream stagnation temperature. The other input parameters for the forebody 
boundary layer Computation (DTDYCA, NTABA, XWLA, TWLA, and DTDYWA) are not used 
and hence are not specified. The input parameters in namel‘*st LISJ.7 used for 
the cowl boundary layer computation .are not entered. 

All input parameters in namelist LISTS retain their default values, thereby 
specifying a case-withro wall mass bleed. The input parameter KDFA is left at. 
its default value, of 0, which specifies an impermeable forebody/centerbody wall. 
The Other input parameters in namelist LISTS used for specification Of the fore- 
body/centerbody mass bleed distribution (NRA, XSA, XEA, and R0VA) do not have 
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to be entered. Since the cowl boundary layer computation is not invoked, the. _ 
•35S bleed input parameters used for the cowl bounda ry lay er need not.be 
entered. 

All input parameters in namelist LIST9 retain their default values except 
for XTl. Consequently, the turbulence model parameters APLUS, XKAPPA, XALPHA, 
PRT, TFACTR, AKLM, and ALAM are left at their default values. The default 
transition model, denoted by ITRANM(1)=5, is retrained. This transition model 
requires that the input parameter array XTl be entered. For this sample case, 
laminar flow has been Specified as existing over the entire length of the fore- 
body. This is accomplished by specifying the onset transition location for 
each circumferential station as being equal to 20.0 ft, hence XT1(1 ,1)-16*20.0. 
The other input parameters in namelist LIST9 do not have to be entered. 

All convergence tolerances and iteration limits retain their default values 
in namelist LISTIO. 

No debug output is to be printed, hence all input parameters in namelist 
LISTll retain their default values. 

Selected portions of the computed output for Sample Case No. 2 are presented 
in Figure 11. 
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RLTLIICIICC VXSCOam» 0.3S0000C«0«fLar-SEC/rT«*3) aerERENCC TENREIIATURr* P.4«2000E^03(DEC fO 




e a G 

5 S 8 

6 2 c ;; 

S O b. 

K» o *n *-* 

M n e c* 

a ct ♦ o 

*<1' « M ♦ 

O f« • M 

» 

O O lA rfl 


e o M 

e at ;» M 

O MX 

« • X a H 

H *« U w iW 

« <« m w 

. 5 “ a s 

as a 8 


»40MM«040«9M« 

kxxxxxxxa 

aoooooooo 


»0 09000 o o 
»e ooooooo 

i)WMUU<«IMWM 

» r> ooooo o o 
io««oooon#t 
» c« o o o r> o 
r <i*oom««e oo 

ro/nmAAmciM 

»o««ooo«o 


g . a : o*. 

X X M • 

O O O X »• 

X *4 M 3 H 

X H K bt »« 

s I H i 

I .- i 8 S 


SSS 88 SaSS;{!;iS 8 S 



Figure IT. Continued. 
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Figure 11. Continued. 
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Figure 11. Continued. 
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4. SAMPLE CASE NO. 3 


Sample Case No. 3 1s concerned with the computation of the internal super- 
sonic core flow for. the Boeing Mach 3.5 mixed-compression inlet described in 
Reference (9). The internal supersonic flow field is computed for the off-design 
conditions of a free-stream Mach number of 2.5, a centarbody forward translation 
of 0.855 ft, and an angle of incidence of 1.0 degree. 

The data deck for Sample Case No. 3 is presented in Figure 12. Again, the 
first card of the data deck is the title card. All input parameters in namelist 
LISTl retain tlieir default values except for KCALL, XEND, KBLAY, and KTRANS. 

Since only the internal flow field is to be computed, KCALL(l) = 0 is specified 
to bypass the forebody flow field computation. The integration termination 
point for the internal flow computation is x=5.05 ft, hence XEND(2)=5.05 is 
entered. For this sample case, neither the boundary layer computation is to be 
performed nor are bleed effects to be accounted for in the computation of the 
supersonic flow; hence, KBLAY>^0 and KTRANS=0 are specified in.namelist LISTl. 

All input parameters in namelist LIST2 retain their default values except for 
MFS, XI, and KSUPER. Since the free-stream Mach number for this sample case is 2.5; 
MFS=2.5 is specified in namelist LIST2. With the prescribed. centerbody translation 
of 0.855 ft, the axial location of the cowl lip, and thereby of the initial -value 
plane, is x=3,715 ft. Note that the origin of the coordinate system remains at 
the forebody rip, hence a forward centerbody translation corresponds to a rearward 
cowl translation. Consequently, XI=3.715 is specified in namelist LIST2. For 
this execution, the results of the approximate analysis described in Section II 
are used for the initial data since the forebody is conical. Hence, the initial- 
value plane flow property field is internally specified by entering KSUPER=1. 

All input parameters in namelist LISTS, which specifies the supersonic 
flow computational mesh, are. retained at their default values. 

All parameters in namelist LIST4, which specifies the thermodynamic and 
molecular transport property models, retain their default values. 

The contours of the centerbody and the cowl are specified in namelist 
LISTS. For this inlet, the centerbody contour and the cowl contour are described 
by equation (3) applied to a number of intervals on each contour. Consequently, 
KBASE retains its default value of 0. The number of axial stations used for 
the geometry description of the centerbody and the cowl are 11 and 14, respec- 
tively; hence NCENT-11 and NC0WL*14. Thus, the number of intervals on the . 
centerbody and the cowl are 10 and 13, respectively. Since equation (3) is - 
used to specify the body radius for all intervals, KDCENT(I)=1 (1=1 to 11) and 
KDC0WL(I)=^1 (1=1 to 14). Note that except for XCENT(Jl) and XC0WLC14), the 
specifications for the other input parameters for the last station on either 
the centerbody or the cowl are immaterial. That is, only the XCENT and XC0WL 
arrays must be.speci fled. for NCENT and NC0WL elements, respectively. The other 
arrays, require that only (NCENT-1) or (NC0WL-1) elements be entered. The 11 
elements of the XCENT array, are thus entered. Although not used in the compu- 
tation, the 11 elements of the RCENT array are entered, each element being 0.0. . 

The arrays ACENT, BCENT, CCENT, and DCENT are -then entered. For each of these 
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arrays, the last (11th) element is arbitrarily specified as 0.0. In a like 
manner, the arrays XC0WL, RC0WL, AC0WL, BC0WL, CC0WL, and DC0WL are entered. 
For this case, the oenterbody has been translated forward by 0.855 ft, or, 
alternati-vely, the cowl has been translated. backward by 0.855 ft. Hence, 
DXTRAN=0.855 is. specified in nameliit LIST5. 


Since in this sample case the boundary layer computation is not invoked, 
the. input parameters for namelists LIST6 and LIST7 need not be entered. 

Mass transfer effects are not Considered for this execution, thereby all 
input parameters in namelist LISTS are retained at their default values. Thereb y , 
impermeable forebody/centerbody and cowl contours are specified. 

The input parameters for namelist LIST9, which is used for specification 
of the turbulence model, need not be entered for this sample case. 

All convergence tolerances and iteration limits specified in namelist-.- 

LISTlO retain their default values. 

No debug output is to be printed, thereby all input parameters in namelist 
LISTll retain their default values. 

Selected portions of the computer printout for Sample Case No. 3 are 
presented in Figure 13. 
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Figure 13. Selected output for Sample Case Ne. 
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Figure 13. Continued. 
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Figure 13. Cor>tinued. 
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Figure 13. Continued. 
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Figure 13. Concluded. 


5. SAMPLE CASE NO. 4 


This sample case is concerned with the computation of both the supersonic 
core flow and forebOdy/centerbOdy and cowl boundary layer flows for the Mach 
3.5 mixed-compression inlet considered in Sample Case No. 3. The internal flow 
field is again computed for the off-design conditions of a free-stream Mach 
number of 2.5, a centerbody forward translation of 0.855 ft, and an angle of 
attack of 1.0 degree. The computation for this sample case includes the effects of 
distributed mass bleed. Independent turbulent transition models are specified 
for the forebody/centerbody and cowl boundary layer calculations. 

The data deck for Sample Case No. 4 is presented in Figure 14. The first 
card of the data deck is again the title card. Namelist LISTl for this sample 
case is identical to thatin Sample Case No. 3 except for the input parameter 
KBLAY. Since the boundary layer computation is to be invoked for this sample 
case, KBLAY retains its default value of 1. Mass transfer effects will not 
be considered in the computation of the supersonic core flow, hence KTRANS 
is retained at the. value of 0. 

Namelists LIST2, LIST3, LIST4, and LI3T5 for this sample are identical to 
those in Sample Case No. 3. 

All input parameters in namelist LIST6 retain their default values except 
for the coordinate stretching parameter FAC.T0R. Hence, the laminar forebody/ 
centerbody boundary layer initial data will be generated internally using the 
Adams algorithm (KBLIDA*=1 and KTURB=0). Since CLENGH=4.0, the length of the 
conical forebody section is specified as 4.0 ft. The Adams algorithm mesh 
parameters ADY and ARATI0 retain their default values of 0.010 and 1.0630. The 
laminar cowl boundary layer initial data are internally generated (KBLIDB=1 anc 

KTURBB=0),. Since MBLAY, NA, and NB retain their default values. Of 15, 20, and 

20 respectively, both the forebody/centerbody and cowl boundary layer computa- 
tions employ 15 circumferential stations and 20 radial stations.. For this 
sample case, the boundary layer normal coordinate stretching factor is specified 
as 1.175 for both boundary layers i. hence, both elements of FACT0R are specified 
as 1.175. 


All input parameters in namelist LIST7 are retained at their default values 
except for TG0NSTA and TC0NSTB. Consequently, constant wall temperature boundary 
conditions are specified for both the forebody/centerbody and cowl boundary 
layer computations (KTYPE-.1, KWLTA^l, and KWLTB=^1). TC0NSTA and TC0NSTB are 
both specified as 876.698, which is the free-stream stagnation temperature of 
876.698 R. 


Distributed mass bleed is specified on both the forebody/centerbody contour 
and the cowl contour for this computation, hence KDFA=1. and KDFB=1 are entered 
in namelist LISTS.. Two bleed zones On the forebody/centerbody contour and one 
bleed zone on the cowl contour are employed, thus NRA=2 and NRB*1. The bleed 
zones for the forebody/centerbody contour extend from x=4.2 ft to x=4.4 ft 

first zone, and from x*4.9 ft to x=5.1 ft for the second zone. Hence, 
XSA(1)=4.2, XEA(1)=4.4, XSA(2)=4.9, and XEA(2)=5.1 are entered in namelist 
LISTS. For the COwl , the bleed zone extends from x=3.75 ft to x=3.F'' ft; hence. 
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XSB(1)=‘3.75 and XEB(1)=3.85. are specified. The bleed mass flux for each zone 
is 0.000001 (slug/ft^.sec). - Consequently, R0VA(1), R0VA(2), and R0VB(1) are 
each entered as 0.000001, 

The turbulence model parameters are specified in namelist LIST9. For 
this computation, transition models No. 5 and No. 3 will be used for the fore- 
body/centerbcdy and cowl boundary layer computations, respectively. Hence, 
ITRANM(1)=5 and ITRANM(2)=3 are entered. Both models require specification 
the transition onset location array XTl. Although a different transition on^ 
set distance for each circumferential Station can be specified, a constant 
distance is used in the present computation. For the forebody/ centerbody, a 
transition onset location of x=4.0 ft is arbitrarily selected; whereas for 
the cowl, a transition onset location of x=0.3 ft is selected. Consequently, 
XT1(1 ,1)=16*4.0 and XT1(1 ,2)-16*0.3 are entered. The transition model used for 
the cowl boundary layer computation requires that the XT2 array be entered to 
denote the onset location of fully turbulent flow. For the present case, all 
circumferential stations on the cowl have x-1.8 ft selected as this fully tur- 
bulent onset location; hence, XT2(1 ,2)=16*1.8 is entered in namelist LIST9. 

All convergence tolerances and iteration limits retain their default values 
in namelist LISTIO. 

No debug output is to be printed for this sample case,_^o all input param- 
eters in namelist LISTll retain their default values. 

Selected portions of the computer output for Sample Case No. 4 are presented 
in Figure 15. 


ORIGINAL PAG!* 13 

OF pnoPv 





II 

\n a 

vt tij 

rt sZ 

H 

btr ^ 
II 


w in 

3 -H 

2 

u.' • 

X rn 

It 

o X 
II 

^ *Jf> 

'T • * 

^ II cs 
•<•>3 2!: 
^ U fj3 
.2 vi: ca 

^ r>i hi 
^ ^ 5 h 
/ i '0 »/l 

X -I M *-l 

.2 4 .4 

a 41 


in 


<rk 


00 

CM 


^ • 

^ tO 

■» 

It ' 


o\ 

rn 


in 

CO 

n 

r- 

• 


m 

Si 


in 

as 

M 


tn 


• tn 


o 

<• 

n 

0 

• 

1 


^ in 

r- 

CM nr 

<• •>. 

JSi ^ C> m 

'X3 CM in 


m 

yc • vn 

Q <• 

• 

a ^ 

iC Si ^ 

in as 

^ • hO 

^ *tt» 

^ m 


It • 
sr 


£-• CM 


r- 

CD 

rn 


o 


m 

o 




in 

o 

o 

o 


as 

St> 

Mn 

CM 

2\ 

m 

0 
• 

1 


o 


o % 
• m 

ca • 


o\ 
o\ 

<* %. 

r^ o 

f^ • 

^ o 


o ^ 


Si tjtZ 

z • 

u: 4* 


h) 

hi 

C OL 

u 

CM 

O 

• -12' ■ 

o % 

hi 


in lo 

.iX o 

in 

«-4 

11 :x 

• 

o 

• 

r*. 

% 4» 

p'* 



^ ~ 

II <* 

i4 as 
^ r- 

O GD 
U 

2S 1^ 


«‘CM 


Si 

rn 

as 

9 

o ^ 

^ ^ o 

II C O o • 

H • n • o 
3 o o . 
-e '*3 *. 

U II II II ^ 
•*a /2 ^ ^ r*» • 

il w w 

H &< H H in 
^0 s} » «: r? ;N 
'-I -H it! aI jn 
2 U «4 O ‘O 
^a a X X <x • 


m 

rn 

n 

r- 


rn 

m 

o 

1^ 


^ r- 

m 00 

m w 

O' »n 

C vO 

CM tH 

O • 

JN 

O • 

• in 

c^ 

o as 

• at 

o ^ 

hi 

^ o 

c • 

• I 


^ o » 
m 

o o • 

0 -^r 
in •. 
r* CM •. 
^ O CM 

• 00 • 

• I fO ^ 
CM 

^ Gt «. 
•N ^ O 
• • 

Si i nr 

r- 

vO % %. 

hi «r 

hi in • 

• Cft hi 

1 CM 

T-l % 

^ m 

• O hi 

o ♦ 


in 

CD 


CM 

Si 

-X> 


•r 

vO 

CD 


JO 

vn 

r^ 

00 


in 

01 


so 

m 

a 


CD 

o 

as 


in 

o 

o 


in 

00 

‘Z> 



^ <-f 

hi 

s> <M 

as 

o ^ m 

o 

O o CM 

•■ • 

o . ^ 

1 

. O O 00 

o 00 

• 

»n ^ o 

hi 

o • 

o 

• 00 1 

T-l 

4* 

• 

^ CM 4* •» 

t 

'O ^ 00 4* 
*jTi *at i4* 1 

*. 

Ot at »-i 

hi 

Ot O h- ir-l 

<«-# 

in • CM # 

CM 

o> t • I 

• . 

CM -t 

1 

• • *. 

1 CM •. O 


•-h in • 

CM 

^ o o o 


Ot m o 

as 

O' (M o 

«h 

O' o o in 

• 

<t ^ rn m 

1 

'XJ o o un 

o • • m 


'<h 1 r-l 

4 

• CO 

CM 

1 ID 


CM o O 

• 

^ at 4* • 

1 

•h ^ t 

o • o 


O 4* ^ - 


O ^ CM 4* 

vt 

o o hi 

O 

in • o at 

* 

4) O • ot 

1 

o i <T^ 


as 


CM 

•r 

041 

a ^ 

U) ^ % 

41 in 
< 00 
00 (0 • 
O' X hi 
XJ II 


hO 

■ - • 

o 

If 

%c 

f-# 

II 

CM 


II ^ 


2 ^ 
a o 

u • 

0 O 


II 

o 
-• o 

w Jl 

CE n 
J ^ 

u • 

4 I 


^ % so 

O hi -Ol 

• o • 

o I (M 

u: 

II >o II 

^•in rS 
^ o ^ 

W 'JV w 

H at J 
'2 1* * % 
'*J • O Ot 

U rM u • 

^ I < n 


o 

o 

» 

•t* 


ID D 
oc « 
^ o 
<* 

o % 

O 

• CM 
*-• CM 

C^ 

^ • 

o 

• « 

•H a 
II hi 
CO 
•M GO 
w ^ 
:2 

% 

Q in 
U no 

< • 


in « 

^ o 

• 

^ o 
O 

• ^ 

a 

%<M 
•-< t^ 

o o 

o • 

in 

r- i> 
o 

O 

• • 


I CM 

O 

^ o 

in in 

h) 

n :> 

n rM 
rM • 
Ot 

■4* ^ 

0 r- 

• at 

1 at 
ot 
«r 

o o 

• hi 

o • 


II ^ II 
^ O vs 00 
^ • ^ ,-< 
O W o 
J .4 0 
^ ^ o 

O O O iO 
U • <J hi 
.O O :j • 


hi • 

0 t 

1 

U3 • 
no Si 

^ o 

<* I 

o :u 

^ 4* 

-f-t 

• 4* ^ a 
4 * «-i o Z 

in • ui 

_•* ."M O. 41 

o- • 

• «-• ^ ;n 
o I dO 

O 00 
II « 1 • 

^ ab ui o 
vf ot Ob n 
^otf^ ar 
t4 at jt < 
at 4* ^ 
Q 4* r^ H 
u rM • X 
O » JO 


vO 

CD 
II 
£0 
H 
Q V3 
z :z 
u a 
Ul o 
H 
tn 

%. 

f-t 00 

• at 

^ vO 

# • 

CM >0 
If r- 
00 
II 


II w 

00 rn 
OS :a 
3. X 


CM Si 
II r- 
< • 
cr hi 
ss II 


H 

X 

Q. 

Z « 

UI o 

41 • 

4 * 
SO 

0 ^n 

1 ^ 
ui II 
o ^ 

• r-t 

II ^ 


II £ 
SC O 
U* X 

a 


ac H 
a ua 
rs ii 
u o 
< u 

lu H 
m i-« 

HI t-| 

a .4 
<4 vll 


• • 

^ in 
II 
< 

U. 4* 
*2 • 
iC 4» 
II 

OD ^ 
H r-t 
cO w 
hi ^ 
i4 ul 
4! X 


^ Eh 
2) X 
> 
o 

DC *4 <4 
^in *43 

to n 41 

0 

1 ^ CD 
U3 w • 

^4 « *4 

• 2 » 3 
vh < ^ 
ff rx «H o3 
CM H II 4L 
II M ^ 

mo 

iHT at - • «-h 

^ rh M 

< to w ^ 

> M CM M 

•D a N - j 
a: 41 X 41 


?4 

'41 


CO 

M 

«4 

41 


128 


Figure 14. Data deck for Sample Case No. 
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Figure 15. Selected output for Sample Case No. 
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Figure 15. Continued. 
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Figure 15. Continued. 
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Figure 15. Continued. 
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Figure T5. Continued 
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Figure 15. Continued. 
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Figure 15 . Continued. 







Figure 15. Continued. 
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